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FOREWORD 


This  interim  technical  report  was  prepared  by  Hughes  Research 
Laboratories  under  Contract  No.  N00014-77-C-0593.  It  describes  work  per- 
formed during  the  14-month  period  from  15  July  1977  to  30  September  1978. 
The  program  manager  is  C.R.  Giuliano.  The  principal  investigator  is 
R.C.  Lind,  who  assumed  this  role  early  in  th'  program  from  V.  Wang. 
Principal  contributors  to  the  program  are  T.R.  O'Meara,  who  has  been 
primarily  involved  in  the  systems  aspects  of  nonlinear  phase  conjuga- 
tion, R.K.  Jain,  who  contributed  to  the  acquisition  and  analysis  of 
detailed  experimental  data;  and  R.W.  Hellwarth,  a consultant  on  the 
program;  and  S.M.  Wandzura,  both  of  whom  made  key  contributions  to  the 
theoretical  understanding  of  nonlinear  phase  conjugation.  In  addition 
to  R.W.  Hellwarth,  A.  Yariv  also  served  as  a consultant  on  this  program. 

The  authors  wish  to  acknowledge  technical  criticism  and  advisory 
support  by  R.L.  Abrams,  manager  of  the  Optical  Physics  Department,  and 
to  thank  T.E.  Horne,  who  assisted  in  the  many  detailed  experiments 
performed  on  this  program.  We  also  wish  to  thank  M.B.  White  of  ONR, 
the  technical  monitor,  and  Col.  M.  O'Neill  of  DARPA  for  their  interest 
in  this  work. 
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SECTION  1 


INTRODUCTION 

This  report  covers  the  first  phase  of  a continuing  program  designed 
to  explore  a recently  recognized  property  of  certain  nonlinear  optical 
interactions  of  generating  conjugate  wavefronts  that  can  be  used  to  cor- 
rect optical  distortions  in  laser  systems.  These  distortions  include 
optical  train  aberrations,  laser  medium  distortions,  and  atmospheric 
propagation  aberrations.  The  program  was  divided  into  three  basic 
areas  that  bridge  the  gaj.  between  a preliminary  exploration  of  the 
applicable  nonlinear  processes  and  realization  of  their  potential  useful- 
ness to  laser  systems.  The  areas  were  (1)  to  measure  quantitatively  the 
properties  of  phase  conjugation  from  stimulated  Brillouin  scattering  (SBS) 
at  0.69  pm,  (2)  to  develop  a theoretical  understanding  of  nonlinear 
phase  conjugation,  and  (3)  to  determine  the  applicability  of  nonlinear 
phase  conju  gation  to  various  systems  of  interest  to  DARPA.  We  made 
significant  advances  in  each  of  these  areas  during  this  program;  these 
accomplishments  are  described  in  this  report. 

A.  BASIC  CONCEPT 

Nonlinear  optical  phase  conjugation  represents  a new  class  of 
optical  interactions  having  the  potential  for  providing  novel  solutions 
to  problems  in  fields  such  as  adaptive  optics  for  high-energy  lasers 
(HELs) , fusion  laser  systems,  imaging  systems,  and  laser  communication 
systems.  Traditionally,  nonlinear  optical  interactions  have  been  ex- 
ploited to  yield  sum  and  difference  frequency  generation,  parametric 
amplification  and  oscillation,  stimulated  scattering,  nonlinear  spec- 
troscopy, etc.  In  the  present  context,  the  desired  output  wave  from  a 
nonlinear  medium  possesses  the  unique  property  of  being  the  spatial 
complex  conjugate,  or  phase  conjugate,  of  the  input  wave. 
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To  introduce  the  concept  of  phase  conjugation  by  nonlinear  optical 
techniques,  consider  the  existing  adaptive-optics  techniques  for  reduc- 
ing the  effects  of  atn’osph  ric  turbulence  and  optical  train  distortions 
on  laser  beams  (coherent  optical  adaptive  techniques,  COAT,  for  example). 
COAT  systems  adjust  the  phase  front  of  the  transmitted  beam  to  compensate 
for  the  phase  distortion  iut.rodu  ed  by  the  optical  medium.  These  sys- 
tems sense  the  distortion  bv  monitoring  backscattcr  from  the  target  and 
form  a phase  front  that  is  a phase  conjugate  of  the  distortion  to  be 
corrected.  The  correction  is  iccomplished  by  using  some  form  of  dis- 
crete, multichannel,  phase-front  corrector  such  as  a deformable  mirror 
driven  by  electronic  servos.  In  contrast,  the  method  of  nonlinear 
optical  phase  conjugation  directly  generates  the  spatial  phase  conju- 
gate of  a distorted  wavefront.  This  phase  conjugate  can  be  transmitted 
through  the  original  distorting  optical  path  to  form  a corrected  beam. 

The  nonlinear  interactions  automatically  perform  a real-time  phase-front 
correction  with  high  spatial  frequency  capability  over  the  entire  cross 
section  of  the  beam  without  any  external  wavefront  sensing  or  electronic 
control s. 

The  basic  system  concept  of  nonlinear  phase  conjugation  for  laser 
beam  correction  through  the  atmosphere  is  shown  schematically  in  Figure  1. 
The  process  involves  the  following  steps: 

• The  first  transmitted  pulse  propagates  to  th<  target 
(assumed  for  this  example  to  be  a single  unresolved 
glint) . The  purpose  of  this  first  pulse  is  to  illumi- 
nate the  glint  target,  which  then  serves  as  a test 
source  at  the  wavelength  of  interest. 

• The  light  reflected  from  the  glint  propagates  through 
the  distorting  medium  towards  the  transceiver,  arriv- 
ing as  an  aberrated  wavefront. 

• The  phase  conjugate  of  this  distorted  re'  irn  wave  is 
generated  by  a nonlinear  optical  conjugator. 

• After  coherent  amplification  to  a desired  power  level, 
the  phase-corrected  pulse  retraverses  the  distorting 
atmosphere,  which  now  restores  its  phase  coherence  so 
that  the  entire  beam  is  focused. 


10 


POWER 


8281-7 


— 0 ILLUMINATION  PULSb 
GLINT  RETURN 

V CORRECTED  HIGH  ENERGY 
PULSE 


Figure  1.  Basic  phase-conjugation  scenario. 
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The  lirst  demon si  rat  ion  of  the  removal  ol  aberrations  in  an  optical 
train  by  nonlinear  phase  conjugation  was  reported  by  Zel'dovich  et  al. 
(JKTP  Lett.  Vo  1 . 15.  pp.  109-113,  1972)  and  Nosach  et  al.  (JK'IT  Lett., 
Vol.  lb,  pp.  435-4)8,  1972).  In  those  experiments,  the  nonlinear 
interaction  used  was  SliS  in  t'S,;.  In  that  case,  the  Brillouin- 
backscat t ered  wave  is  the  complex  conjugate  ol  the  input  signal. 
Recently,  these  experiments  have  been  verified  and  o*  Hied  at  URL 
(see  Appendix  A).  This  initial  work  led  to  the  program  described  in 
this  interim  report. 

B.  SUMMARY  OK  PRINCIPAL  ACCOMPLISHMENTS 

The  main  accomplishments  of  the  first  phase  of  the  program  are 
listed  below: 

• Demonstrated  complete  aberration  correction  by  SBS  for 
beams  aberrated  to  35X  the  diffract i n l'aic. 

• Performed  detailed  photographic  measurements  to 
establish  spatial  distribution  of  beam  proliles  to 
intensities  down  to  MO  of  peak  intensity. 

• Established  experimentally  that  the  fraction  of  nonconju- 
gate return  when  using  SBS  for  correction  is  below 
measurement  Limits. 

• Developed  systems  applications  exploiting  novel  varia- 
tions of  nonlinear  phase  conjugation  techniques  in  the 
area  of  high-power  oscillator  compensation,  master- 
oscillator  power-amplifier  train  compensation,  and  for 
systems  which  require  correction  for  atmospheric  tur- 
bulence with  very  high  spatial  frequency  content. 

• Furthered  the  theoretical  understanding  of  phase 
conjugation  by  SBS  in  waveguides  and  In  uonwaveguide 
( f ree- space)  geomet r ies . 

This  report  is  divided  into  four  sections  followed  by  five 
appendices.  Section  2 describes  the  experimental  program  on  SBS  phase 
conjugation  at  0.69  pm.  Section  3 discusses  theoretical  considerations 
relevant  to  SBS  phase  conjugation  and  describes  phase  conjugation  by 
four-wave  and  three-wave  mixing  techniques.  Section  4 describes  poten- 
tial applications  using  the  four-wave  mixing  phase-conjugatior  scheme. 
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Of  KKlMKNi'Al.  STUUIKS  OK  SHS-I'IIASK  t'ON.HHIAT ION  AT  O.tOI  tun 

Tho  primary  oh|oot  iv«>  ol  t ho  no  oxpot Imonl  a wan  to  ilo t oimlno  t ho 
ampl  i t into  ot  tho  SI'S  haoksoat  1 01  oil  wavo  anil  tin*  aoonraov  ol  oorroot  Ion 
lot  I ho  SI'S  Wiivo  an  a Innot  Ion  ol  paramot ora  anoli  aa  tin*  goomot  i v ol 
tlio  lilt  oraot  Ion,  ilogt  oo  ol  ahorrat  Ion,  ami  powor  ahovo  t hronliolil , 

Wavogn  IiIoh  ot  vnrlonn  longtha  anti  "t  roo  apaoo,"  or  nonwavogn  Itloil , 
goomotrlon  woio  at  ml  toil.  In  addition,  imago  i ooonat  i not  Ion  lining 
pliaao  ooii  |ngat  oil  SIIS  wan  ilomoiiwt  rat  ml . Tho  donor Ipt ion  ot  thoso  ox- 
porlmonta  la  pi  ooodoil  hv  a hi  lot  ii  laonaslon  ot  phaao  oonliigat  ion  hv  HUS. 

A.  I'll  ASK  tmilHIATION  I'Y  SltS 

An  umloiHt  ami  lug  ol  tho  phvaloa  ol  tho  SIIS  protons  la  miloil  hv 
totalling  that  light  la  aoattotod  I tom  an  aoonat  io  wavo  in  a maitnor 
lilontloal  to  tho  aoatloiing  t tom  a ilittiaot  Ion  grating.  Tho  mot  ion 
ol  tho  aonml  wavoa  toaiilta  in  a Dopplor  aliitt  ol  tho  station'll  light. 

In  SUS,  a at  tong  otootiio  tlolil  ‘an  ho  ptoilnooil  hv  tho  paaaago  ot  an 
intouno  light  hoam.  Thtough  oloot tost  riot  ion,  thin  toaiilta  In  portoillo 
ohangoa  in  tho  ilonaltv  ol  tho  moilliiin  ami  thoroloro  in  tho  moillnm  In- 
ilox.  whloh  gonotaloa  a ttavollng  aoonat  lo  wavo  whitli,  In  turn,  aoattoia 
(tolloota)  aomo  ol  tho  input  opi ioal  hoam. 

Klgmo  ptoaonta  a almplltloil  vlow  ol  thin  protons,  Tho  phaso- 

a hot  tat  oil  hoam  la  topi  oaont  oil  hv  a wavo  t tout  with  a almplo  atop.  An 

otillnntv  lot  loot  ion  toaiilta  In  tho  ilonhlv  ahoitatoil  w.ivotront  aliown, 

Aa  a toanlt  ot  tho  ShS  intotaot  ion,  tho  ahori'al  oil  pump  wavo  otoatoa 

ahoitatoil  aoonat  lo  wavoa,  whloh  aot  aa  a moving  illoloottlo  tol  lootoi  to 

vlolit  tho  ooningalo  aoattotoil  wavo.  It  tho  gonorat  ion  ol  aonml  in  SI'S 

N la  vlowoil  aa  a ool  I ialonal  ptoooaa  In  whloh  a phot  on  (k  ,,u',t)  nililorgooa 

paramot  t lo  tiooay  Into  a phonon  (k  ,,u  ) plna  a photon  (k.,u>  ),  oonaorva- 

a a I I 

tlon  ol  onoi  gv  ami  moitiont  nm  illtoolly  vlolil  (hh  .shown  in  l-'ignto  \) : 
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Figure  2.  Simplified  view  of  the 


SBS  phase-conjugation  process 
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Figure  3.  Generation  of  sound  by  parametric  decay  in  SBS  of  light. 
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The  latter  expression  is  known  as  the  Bragg  condition.  Using  Eq.  1 
and  (i)^  ■ c |kjJ,  we  have 


V 

s 


(3) 


where  n is  the  index  of  refraction  (assumed  to  be  the  same  for  and  (a2) 
and  is  the  speed  of  sound  in  the  medium.  Taking  absolute  magnitudes 
of  Eq.  2,  we  have 

|kj2  - |k2| 2 + IkJ2  - 2|k1|  |k2|  cos  012  . (4) 


For  a collintar  interaction  0 , the  angle  between  waves  1 and  2 is  ir, 

and  cos  0 “ “1.  Hence,  we  obtain 

I kg  I 2 - |k2|2  + |kx|  2 + 2|kJ  |k2|  - (|kj  + | k2 1 ) 2 (5) 

and 


(6) 


By  dividing  Eq.  3 by  Eq.  6 and  rearranging,  we  obtain 


and,  since  n V /c  <<  1 10  ^),  it  follows  that  |k  | » |k  | ■ k and 

we  obtain 
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From  the  appropriate  wave  equations  for  the  photons  and  phonons, 
the  scattered  wave  (k  ,ui. ) and  the  acoustic  wave  (k  ,u)  ) can  both  be 

* 1 S 8 

shown  to  experience  gain.  If  the  gain  is  greater  than  the  loss  in  the 
medium,  a growing  acoustic  and  scattered  light  wave  will  he  produced. 

Ihe  condition  for  highest  gain  is  for  collincar  waves  with  the  scattered 
light  traveling  counter  propagating  to  the  input  beam.  Using  the  con- 
dition for  the  gain  to  exceed  the  loss,  the  threshold  pump  intensity 
(1^)  is  =800  MW/cm  for  CS^*  This  result  is,  as  a first  approximation, 
independent  of  the  wavelength  of  the  input  photons.  The  frequency 
shift  for  these  cases  is  calculated  to  be  Aca/to  - 10~5,  or  ^ 6 GHz,  at 
6943  A for  CS2  and  Aw/w  ^ 1.3  x Iff4,  or  4 GHz,  at  10.6  pm  for  Ge. 

At  the  high  power  densities  necessary  to  reach  threshold  for  SBS, 
care  must  be  taken  to  avoid  self-focusing.  The  threshold  for  self- 
f ocusing  is  dependent  on  total  power,  varies  as  a function  of  wave- 
length squared,  and  has  a value  of  % 20  kW  for  CS2  at  ruby  wavelengths. 
Ihese  basic  parameters  are  relevant  to  the  experiments  to  be  described 
below. 

B.  EXPERIMENTAL  APPARATUS 

A schematic  of  the  basic  experimental  configuration  is  shown  in 
figure  4,  A ruby  oscillator/ amplif ier  combination  generates  an  output 
beam,  which  is  focused  either  into  a multimode  waveguide  of  a pre- 
scribed length  filled  with  CS^  or  into  a large  cell  containing  CS2* 

The  use  oi  a waveguide  resuits  from  our  early  experiments  (see 
Appendix  A)  verifying  the  original  results  of  Zeldovieh  et  a!.2  One 
of  the  issues  addressed  on  this  program  has  been  to  determine  the  ef- 
fects of  waveguide  and  nonwaveguide  geometries  on  the  phase  conjugation 
process  (see  below).  The  SBS  backward-generated  wave  is  split  off  by  a 
beam  splitter  and  focused  at  the  film  plane  with  a 1-m-f ocal-length 
lens  for  observation.  A detailed  diagram  of  the  apparatus  is  given 
in  figure  S,  Ihe  ruby  laser,  passively  Q-switched  by  a dye  cell. 
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Figure  4.  Experimental  eon!  igurat  ion  . 
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Figure  5.  Detailed  diagram  of  0.69-pm  experimental  apparatus. 


provides  a TEM^  single-longitudinal-mode  output  pulse  of  17-nsec 
duration.  This  is  amplified  and  controlled  by  an  adjustable  attenuator 
to  deliver  from  10  to  100  mJ  to  the  SBS  cel] . 

An  aberration  can  be  purposely  introduced  into  the  beam  at  the 
position  shown  in  Figure  5 to  test  the  corrective  properties  of  SBS 
backscatter.  The  aberrator  is  a microscope  slide  etched  in  HF.  The 
original  beam  quality  is  monitored  by  removing  the  lens-CS^  cell  combi- 
nation, placing  a partially  reflecting  mirror  in  that  position,  and 
reflecting  a portion  of  the  light  back  into  the  1-m-focal-length  camera 
through  the  beamsplitter.  The  degree  of  aberration  introduced  is  then 
measured  by  placing  the  aberrator  in  front  of  the  reference  reflector, 
which  causes  the  light  to  double-pass  the  aberrator.  The  SBS  correc- 
tion process  is  accomplished  by  passing  the  light  through  the  amplifier 
and  attenuator,  through  the  aberrator  shown  in  position  in  Figure  5, 
and  into  the  focusing  lens  and  CS^  cell.  For  the  waveguide  experiments 
the  lens  is  used  to  focus  the  beam  to  a waist  just  in  front  of  the  cell 
entrance  window,  giving  a divergent  entrance  beam.  Since  the  index 
of  CS2  is  higher  than  the  glass,  this  diverging  beam  is  completely 
internally  reflected  in  the  capillary.  The  waveguide  SBS  cell  serves 
to  increase  the  interaction  length  and  percentage  of  light  back- 
scattered  without  increasing  the  power  density  to  the  point  where  laser 
induced  breakdown  of  CS^  can  occur  and  without  increasing  the  total 
power  input  to  the  point  of  self-focusing.  The  slightly  red-shifted 
SBS  backscatter  retraces  the  path  of  the  original  beam  through  the 
multiple  bounces  in  the  waveguide  and  through  the  aberrator,  where  a 
portion  is  sampled  by  the  beamsplitter  and  1-ra  camera.  For  the  free- 
space  (i.e.,  nonwaveguide)  experiments,  the  waveguide  is  replaced  by  a 
cell  filled  with  CS^  (a  10-cm  spectrophotometer  cell  is  used).  Various 
focal  length  lenses  are  used  to  focus  the  beam  into  this  cell,  taking 
care  to  avoid  self -focusing. 


(Experiments  were  performed  with  a single  aberrator  (as  shown)  and 
also  with  two  separate  aborrators  placed  at  several  positions  between 
the  attenuator  and  t ho  cell.  The  later  set-up  was  used  to  simulate 
a distributed  aberration,  such  as  encountered  in  a turbulent  atmosphere. 

The  laser  pulse-shape  and  SBS  backseat tered  pulse  were  measured  on 
every  shot  using  the  fast  photodiodes  shown.  In  addition,  the  total 
amplified  input  laser  energy  was  monitored  on  every  shot.  For  measuring 
the  SHS  frequency  shift  beam,  a Fahry-Perot  etalon  was  used  in  the 
position  shown.  For  these  measurements,  a reference  beam  is  also  pro- 
vided for  comparison.  This  was  achieved  by  inserting  a beam  splitter 
in  front  of  the  57-mm  lens  and  directing  the  beam  into  the  Fnbry-Perot. 
Far-field  intensity  functions,  or  divergences  of  the  beam,  are  measured 
at  the  local  point  of  the  1-m  lens  by  a photographic  technique  using 
high-contrast  film  in  a series  of  exposures  at  differing  attenuations. 

In  this  way,  the  actual  profile  of  the  original  aberrated  and  corrected 
beams  can  be  determined  independently  of  film  linearity.  Figure  b 
illustrates  this  technique.  An  arbitrary  beam  profile  is  shown  for 
various  intensity  levels  on  the  film.  The  film  threshold  intensity  is 
indicated  on  the  figure.  The  procedure  followed  consists  of  exposing 
the  film  with  intensity  I and  measuring  the  diameter  at  the  film 
threshold  D.  The  beam  is  then  attenuated  a known  amount  and  the  result- 
ing smaller  diameter  is  measured.  This  procedure  is  repeated  several 
times  (typically  10  to  15)  until  the  intensity  is  reduced  so  that  the 
spot  is  just  resolvable  on  the  film.  In  this  manner,  the  entire  beam 
profile  is  traced  out  independently  of  Him  linearity.  Using  this 
technique,  we  have  measured  the  profile  far  out  into  the  wings  of  the 
distribution  down  to  intensities  of  10  of  the  peak  intensity  (see 
below).  Measurements  were  made  under  the  following  conditions:  (1) 
laser  only  (the  lens-waveguide  combination  was  replaced  with  the  100% 
reflector  in  place),  (2)  SBS  backseattered  beam  without  aberrator  in 
place,  and  (3)  SBS  backseattered  beam  with  aberrator  in  place  and 
double  passed. 


21 


6867— 5R2 


Figure  6.  Determination  cf  beam  profile. 
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EXPERIMENTAL  RESULTS 


An  example  ol  the  type  of  photographic  data  generated  is  shown  tu 
Figure  7.  The  7(a)  through  7(e)  sequenee  is  aherrator  data  without 
SBS:  (a)  the  uuaherratod  beam  (which  is  the  same  as  the  retro- 

ref  lected  beam),  (b)  the  aberrated  beam  passing  through  an  aherrator 
labeled  No.  2 (which  gives  a beam  ^ 8X  diffraction  limited,  and  (e)  the 
aberrated  beam  passing  through  aberrators  No.  2 and  No.  3 (which  gives 
a beam  'v  35X  diffraction  limited).  Figure  7(d)  has  the  beam  passing 
through  the  aberrators  as  before,  but  now  the  phase-conjugate  SBS  beam 
is  generated,  and  the  beam  retraverses  the  aberrators.  In  this  case, 
the  spot  is  essentially  identical  with  either  aherrator  No.  2,  with 
both  aberrators,  or  without  an  aberralor.  These  data  imply  that  the 
SBS-corrected  beam  has  a fat-field  divergence  nearly  equal  to  that  of 
the  input  laser  beam  or  a bear,,  that  is  aberrated  to  3rjX  diffraction 
limited  c.-.u  be  corrected. 

Figure  8 gives  a qualitative  picture  of  the  distortion  introduced 
by  the  aberrators  employed.  This  figure  compares  a photograph  of  a 
resolution  chart  taken  through  aherrator  No.  3 with  a photograph  taken 
without  the  aherrator  present.  A more  quantitative  character izat ion 
of  the  aberrators  used  in  the  above  experiments  has  also  been  made. 

The  aberrators  were  scamud  mechanically  by  running  a stylus  across 
the  surface  of  each.  An  example  of  such  a scan  tor  aherrator  No.  3 is 
shown  in  Figure  9.  The  peak-to-penk  phase  variation  over  a distance 
corresponding  to  the  ruby  beam  diameter  is  typically  a *2?i  phase  shift. 
If  the  data  from  such  scans  were  digitized,  statistical  information, 
rms  deviations,  etc.  could  be  obtained.  Also,  fast  Fourier  transform 
techniques  will  give  spatial  frequency  data  for  the  aberrators.  Such 
additional  data  might  be  useful  if  a detailed  numerical  description  is 
required. 

Using  the  photographic  data,  two  methods  of  displaying  the  results 
can  be  used.  First,  a semi-log  plot  of  relative  intensity  versus  spot 
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diameter  squared  as  recorded  on  t lie  film  could  be  made.  A straight-1 lue 
fit  will  result  if  the  beam  is  Gaussian  with  a slope  proportional  to  the 
1/e  (In  intensity)  Gaussian  beam  radius.  To  see  this,  consider  the 
following,  begin  with 
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where  l'  Is  the  total  power  in  the  beam,  und  iS  is  the  l/e  radius  (in 
intensity)  of  the  beam.  If  multiple  exposures  are  taken,  we  obtain 
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Assume  that  a particular  total  power  and  beam  radius  6 (common  to 
all  spots)  give  1^,  (the  Intensity  required  to  reach  the  film  exposure 
t hresho Id) : 
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For  each  spot  measured,  we  know  P,  , l’  ...  ; thus, 
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and  we  have 


which  is  the  desired  result.  The  validity  of  this  technique  is  based 
on  the  assumption  that  the  film  possesses  a sharp  exposure  threshold 
and  that  the  diameter  of  a given  spot  is  equal  to  the  beam  diameter  at 
which  the  intensity  equals  the  exposure  threshold.  Semi-log  plots  of  this 
kind  are  given  in  Figure  10.  Figure  10(a)  gives  the  aberrated  beam,  which 
is  a good  approximation  to  a Gaussian  beam  of  about  11  times  the 
divergence  of  the  original  beam.  Similarly,  Figure  10(b)  plots  the 
nearly  diffraction  limited  SBS-corrected  beam.  The  divergences  shown 
are  calculated  from  the  above  formula  using  a straight-line  fit  and  a 
1-m-focal-length  lens.  This  type  of  data  reduction  clearly  indicates  that 
corrected  beams  can  be  achieved. 

In  addition  to  this  data  and  all  those  to  be  shown  below,  many  addi- 
tional runs  were  made.  These  were  performed  using  various  aberrators  in 
various  combinations  and  locations  in  the  beam  path.  In  all  cases, 
essentially  perfect  correction  was  achieved. 

The  analysis  of  these  data,  however,  does  not  determine  the 
intensity  fraction  of  the  SBS  beam  that  is  the  phase  conjugate  of  the 
input  wave  and  the  fraction  that  is  not.  This  can  only  come  from  the 
detailed  photographic  data  coupled  with  power-in-the-bucket  measure- 
ments to  be  described  below.  To  see  why  this  is  important,  consider 
the  example  profiles  drawn  in  Figure  11.  Two  profiles  are  shown: 
the  incident  beam  and  the  SBS-corrected  beam.  For  this  example,  the 
intensity  of  the  SBS  beam  is  10%  of  the  incident  beam,  a value  typical 
of  the  present  experiments.  The  dashed  line  corresponds  to  the  film 
threshold  intensity.  Assume  that  the  SBS-corrected  profile  follows  the 
solid  line  as  shown  to  intensities  well  below  the  film  threshold.  In  this 
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Figure  11.  Example  profiles. 


east',  two  facts  emerge.  First,  the  SBS  profile  has  a shape  Identical 
to  the  incident  profile.  Second,  all  ol  the  energy  in  the  beam  Is 
contained  within  this  profile.  This  implies  that  the  SBS-eorrect ed 
beam  is  the  phase-conjugate  beam  and  that  the  fraction  of  conjugate 
return  is  100%.  Now  assume  that  the  SBS  beam  has  the  profile  given  by 
the  dotted  line,  clearly  the  worst  ease  possible  and  unlikely  to  occur. 
In  this  case,  a large  fraction  of  the  energy  could  be  found  in  the 
region  below  the  dotted  line  (this  would,  of  course,  be  non-conjngate- 
veturn  energy).  The  objective  of  the  detailed  measurements  Is  then 
twofold:  to  measure  the  Intensity  as  far  into  the  wings  as  possible 

and  to  make  at  least  one  pinhole  measurement  of  the  energy  to  confirm 
its  distribution.  This  has  been  done  for  several  experimental  con- 
figurations; the  detailed  results  are  presented  below. 

From  a complete  sequence  of  photographic  data,  Figure  12  shows  a 
normalized  intensity  profile  as  a function  of  beam  radius  for  a CS.}- 
f l i led  waveguide  cell  81  cm  long  x 2.!>  mm  in  diameter.  Data  are  given 
for  the  retroref looted  beam,  SBS  with  no  aberrator,  SBS  with  aberrator 
No.  2,  and  SBS  with  aberrators  No.  2 and  No.  3.  To  plot  the  data  in 
this  way  requires  an  additional  data  reduction  technique.  All  se- 
quences ar“  normalized  to  the  same  point,  which  requires  a common  in- 
tensity value  for  all  sequences.  This  value  is  the  exposure  threshold 
intensity.  The  exposure  threshold  is  determined  by  plotting  intensity 
versus  diameter  squared  and  extrapolating  the  straight-line  fit  to  the 
point  of  zero  diameter.  This  procedure  is  shown  in  Figure  13  for  all 
the  data  runs  shown  in  Figure  12  and  those  to  follow.  Figure  12  shows 

(as  discussed  previously)  that  the  SBS-correct ed  beam  has  the  same  pro- 

-3 

file  as  the  uuaberrated  beam  down  to  'V  ]0  of  the  peak.  Pinhole 
measurements  were  made  for  the  conditions  of  Figure  12  using  the  pinholi 
diameters  shown  end  for  an  additional  larger  diameter  pinhole  that  col- 
lects all  of  the  energy,  (A  description  of  the  apparatus  for  pinhole 
measurements  is  discussed  below.)  By  integiatlng  the  profiles  shown 
and  comparing  these  values  to  the  pinhole  measurements  we  can  show  that 
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fraction  of  nonconjugated  return  produced  is  immeasurably  small.  That 
is,  energy  that  could  contribute  to  a return  that  deviates  from  the  pro- 
file shown  (e.g.,  non-phase-conjugate)  is  zero  to  within  the  experi- 
mental error  (estimated  to  be  'v  15%). 

-4 

Similar  results  (measurements  down  to  ^ 10  of  the  peak)  have 
been  obtained  for  waveguides  shorter  than  81  cm.  Figure  14  gives  such 
data  for  an  11-cm,  a 6-cm,  and  a 3-cm  tube  (the  81-cm  data  is  also 
given).  A comparison  of  these  results  again  indicates  good  correction, 
as  in  the  81-cm  case,  but  with  a slight  trend  at  the  lower  intensities 
toward  larger  divergences.  This  indicates  that  a larger  fraction  of  the 
energy  may  be  in  the  wings  for  the  shorter  tubes  than  for  the  81-cm 
tube.  The  total  SBS-backscattered  energy  for  all  waveguide  lengths  was 
measured  to  be  'v  10%  for  cases  without  an  aberrator  and  with  aberrator 
No.  2 double  passed. 

The  pinhole  measurement  apparatus  is  shown  in  Figure  15.  Using 
the  beamsplitter,  a film  plane  and  cross-hairs  were  located  at  the 
focal  plane  of  the  lens  with  the  pinhole  also  at  the  focal  plane.  A 
photograph  of  the  spot  was  taken  for  each  datum  point.  This  spot  was 
compared  to  the  cross  hairs  to  observe  if  a symmetric  illumination  of 
the  cross  hairs  can  be  achieved.  A symmetric  illumination  implies 
that  the  beam  is  entering  the  center  of  the  pinhole.  These  data  were 
then  recorded. 

Experiments  were  conducted  to  determine  the  effectiveness  of 
phase  conjugation  as  a function  of  power  above  threshold.  The  81-cm 
tube  was  used  to  perform  the  experiments.  The  far-field  beam  diverg- 
ence was  measured  over  input  powers  ranging  from  about  10%  over  the 
SBS  threshold  to  four  times  the  SBS  threshold  The  far-field  beam  di- 
vergence obtained  in  these  measurements  was  independent  of  pump  power, 
indicating  that,  for  a given  waveguide  geometry,  the  degree  of  conju- 
gation is  independent  of  power  over  the  SBS  threshold. 

In  addition  to  the  experiments  described  above,  which  were  con- 
ducted in  a waveguide  geometry,  experiments  were  performed  to  study 
the  possibility  of  phase-conjugation  via  SBS  in  a fret-space  geometry. 
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A 4. 5-in.-f oral-length  lens  was  used  to  torus  the  ruby  laser  beam  in  a 
large  cell  filled  with  CS,,.  The  experimental  arrangement  was  essentially 
the  same  as  that  shown  In  Figure  4,  and  the  same  photographic  techniques 
as  those  used  in  the  waveguide  experiments  were  used  to  analyze  the 
tree-space  data. 

In  contrast  to  the  waveguide  experiments,  the  backward-scattered 
SBS  signals  showed  only  partial  wavefront  correction  and  relatively 
Incomplete  phase  conjugation  when  the  pump  power  was  only  a small 
fraction  above  the  threshold  for  SBS.  However,  the  degree  of  conjuga- 
tion was  observed  to  depend  critically  on  power,  and  every  small  incre- 
ment in  the  pump  power  resulted  in  a significant  improvement  in  the 
aberration  removal  capability  of  the  backscattered  SBS  wave.  In  addi- 
tion, the  degree  of  phase  conjugation  improved  significantly  with  the 
introduction  of  random  aberrations  or  strong  transverse  random  ampli- 
tude fluctuations  in  the  beam  (with  the  scattering  ot  the  pump  laser 
beam  into  higher  order  spatial  modes). 

Figures  lb,  17,  and  18  illustrate  the  strong  dependence  of  the 
degree  of  correction  (phase-conjugat * on)  on  both  the  pump  power  and  on 
the  presence  of  random  amplitude  fluctuations  (i.e.,  aberrations)  in 
the  pump  beam.  The  solid  lines  marked  R34  indicate  the  beam  profile 
of  the  unaberrated  retro-ref lected  beam.  Tbe  data  In  Figure  16  cor- 
respond to  a pump  power  level  (designated  as  i ) that  is  *VS  times  the 

o 

SBS  threshold  power.  At  this  power  level,  the  unaberrated  pump  beam 
results  in  a nearly  phase-conjugate  backward-scattered  SBS  beam  that 
has  v 30  to  40%  of  the  incident  pump  power.  However,  as  implied  by 
the  broader  beam  profile  (R32,  crosses),  this  backward-going  beam  shows 
incomplete  reconstruction,  which  is  evidence  of  Incomplete  phase  con- 
jogat  ion.  introduction  of  the  aberrator  (aberrator  No.  2)  in  the  pump 
path  results  in  almost  complete  collapse  of  the  backward-going  beam 
profile  to  that  obtained  for  retro-reflection,  implying  almost  com- 
plete phase-conjugat ion  for  this  case.  For  this  case,  typical  backward- 
wave  returns  wore  between  S and  10%. 
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Figure  16.  Phase-conjugation  via  free-space  SBS 
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As  seen  in  Figures  17  and  18,  only  slight  increases  in  power 
(10%  and  32%,  respectively)  result  in  significant  improvements  in  the 
reconstructed  beam  (curves  R54  and  R56) . Higher  powers  were  not  useo 
for  fear  of  exceeding  the  damage  threshold  of  the  window  of  the  (JS 
cel). 

The  results  described  above  on  phase  conjugation  in  free  s,  ace 
are  in  complete  agreement  with  the  previously  referenced  paper  by 
Zeldovich  et.  al.  on  the  nature  of  phase  conjugation  via  SBS  and  ar^.  the 
first  quantitative  experimental  results  on  this  behavior.  They  indicate 
the  freedom  from  the  restraint  of  using  a waveguide  cell  in  situations 
where  the  lack  of  beam-pointing  accuracy  might  inhibit  efficient  coupl- 
ing into  a waveguide  cell. 

Figure  19  shows  measurements  made  on  the  frequency  of  the  return 
beam.  A Fabry-Perot  etalon  having  a 3.16  mm  spacing  was  used.  The 
expected  SBS  frequency  shift  of  6 GHz  in  CS^  was  observed. 

Image  reconstruction  via  phase-conjugated  SBS  was  also  demon- 
strated on  the  present  program.  In  these  experiments,  SBS  is  obtained 
by  irradiating  CS^  (see  Figure  20)  in  a waveguide  cell  from  the  Q- 

switched  ruby  laser,  as  in  the  previous  experiments.  The  TEM  -mode 

oo 

ruby  laser  output  passes  through  a transmitting  resolution  chart  to 
acquire  a well-defined  spatial  variation  in  the  transverse  direction 
(Figure  21).  A 200-mm  lens  focuses  this  beam  onto  the  2. 5-mm-diameter , 
81-cm-long  CS2  waveguide  cell,  which  acts  as  a "Brillouin  mirror."  It 
generates  a backward-going  phase-conjugated  wavefront,  which  recon- 
structs the  spatial  profile  of  the  resolution  chart  back  at  the  origi- 
nal plane  of  the  chart  (Figure  21(c)).  By  interposing  a beam-splitter 
between  the  focusing  lens  and  the  resolution  chart,  the  backward- 
going (pseudoscopic)  image  is  reconstructed  at  a second  plane.  With 
the  use  of  appropriate  optics,  this  pseudoscopic  image  plane  is  re- 
imaged on  the  film  of  the  camera. 

When  a severe  aberrator  is  placed  in  the  path  of  the  spatially 
modulated  beam,  it  loses  all  its  spatial  information  in  the  far  field, 
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Figure  19.  Free-space  SBS  frequency-shift 
measurement . 
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Figure  20.  Experimental  arrangement  — image  reconstruction. 
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as  can  be  observed  (Figure  21(b))  by  replacing  the  Brillouin  mirror  with 
a plane  reflector.  However,  when  the  Brillouin  mirror  is  re-introduced 
in  the  experimental  arrangement,  it  phase  conjugates  the  wavefront  of 
the  severely  aberrated  beam  so  that  the  backward-going  beam  undoes  its 
phase  distortion  as  it  retraverses  the  aberrator,  thereby  reconstructing 
the  original  spatial  information  at  the  pseudoscopic  image  plane 
(Figure  21(d)). 

With  an  aberrator  that  results  In  an  approximately  10-times- 
d If  fraction-limited  beam  in  the  far  field,  a spatial  resolution  of 

5 lines/mm  has  been  obtained  with  an  optical  system  having  a limiting 
resolution  of  'V  20  lines/mm.  This  demonstrates  the  potential  of  SBS 
for  image  reconstruction. 

Although  the  degree  of  phase  conjugation  theoretically  obtainable 
from  SBS  may  suffer  from  limitations  (see  below)  not  present  in  other 
techniques  (such  as  four-wave  mixing),  SBS  does  possess  the  advantage 
of  experimental  simplicity.  As  a consequence,  it  may  be  a feasible 
technique  for  phase  conjugation  and  real-time  imaging  applications. 

The  limitations  in  the  spatial  resolution  obtainable  with  SBS  phase 
conjugation  of  severely  aberrated  images,  which  are  difficult  to 
establish  theoretically,  were  elucidated  by  these  experimental  studies. 
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SECTION  3 


THEORETICAL  STUDY  OF  NONLINEAR  PHASE  CONJUGATION  SCHEMES 

A.  SBS  AND  SRS 

A detailed  theoretical  analysis  of  wavefront  reproduction  by  SBS 
and  stimulated  Raman  scattering  (SRS)  in  a waveguide  was  performed  on 
this  program  (see  Appendix  B) . The  solution  obtained  was  for  the 
backward  optical  wave  stimulated  by  a multimode  incident  optical  wave 
in  a waveguide  filled  with  a transparent  nonlinear  medium.  The  basic 
conclusion  from  the  theory  is  that,  under  many  conditions  (e.g., 
number  of  excited  modes,  phase  distribution  between  the  modes,  random 
distribution  of  the  amplitude  of  the  modes) , there  is  one  back- 
scattered  wave  that  has  an  exponential  gain  about  twice  that  of  any 
other  solution  and  that  this  wave  is  essentially  a "phase  conjugate" 
of  the  incident  wave.  The  fraction  of  energy  that  is  in  waves  other 
than  this  phase-conjugate  wave  is  smalls  it  was  found  to  be  no 
greater  than  ^ 7%  for  the  arbitrary  pump  waves  studied.  For  SBS,  it 
was  found  that  the  waveguide  length  should  be  L > SA  1N  1 (where  S - 
area  of  guide,  X * wavelength,  N ■ number  of  excited  modes)  for  the 
phase-conjugate  solution  to  exhibit  the  above  properties.  SRS  was 
also  investigated.  Here  the  length  requirement  is 


, < — 2 

SP.S  'v  N | AX  | 

where  rQ  is  a parameter  approximately  equal  to  'v  0.1,  and  |AA|  = 

I ^stimulated  ^ input  ^ * 

The  theory  also  suggests  that  SBS  would  yield  conjugate  waves  in 
free  space  and  in  multimode  waveguides.  This  was  observed  in  our 
experiments. 
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In  addition  to  understanding  the  basic  theory  of  SBS,  the  general 

question  of  compensation  of  optical  path  difference  (OPD)  errors  rather 

than  optical  phase  errors  (this  distinction  only  appears  in  a frequency 

shifting  device)  was  investigated.  Although  SBS  is  likely  to  correct 

phase  error  rather  than  path  error,  the  potential  existence  of  a path- 

length-error-correcting  device  might  be  important  in  some  applications 

requiring  larger  frequency  shifts.  The  error  A<f>  introduced  by  phase 

correction  with  a frequency  shift  is  small  if  Au)/u)  < 1 (neglecting 

dispersion).  This  criterion  is  easily  met  by  SBS,  where  Aw/w  « 10 

(A<p  is  a typical  relative  phase  error  across  the  aperture).  For  SRS, 

_2 

where  (Aui/uj)  = 10  , apparently  perfect  correction  has  been  reported 

in  the  Soviet  literature  for  phase  aberrations  of  - 2ir  (an  extreme 
aberration  in  terms  of  reduction  in  on-axis  intensity).  When  frequency 
shifts  of  10%  or  greater  are  desired,  phase  aberrations  of  a few  wave- 
lengths would  make  the  question  of  OPD  correction  important.  It  ap- 
pears that  nonlinear  phase  conjugating  processes  that  are  presently  con 
ceived  are  not  OPD-correcting  devices.  II  such  an  OPD  device  were 
possible,  it  would  exhibit  several  anomalous  characteristics  and  would 
not  be  capable  of  correcting  all  forms  of  phase  aberration.  This  is 
traceable  to  the  periodic  nature  of  any  phase  variable.  The  anomalous 
characteristics  Include  (1)  nonuniqueness  of  output  for  a given  input; 
(2)  lack  of  linear  superposition,  giving  rise  to  image  dependent  dis- 
tortions; and  (3)  (spatially)  discontinuous  outputs  for  certain  con- 
tinuous inputs. 

B . FOUR-WAVE  MIXING 

A detailed  theoretical  development  of  a second  phase-conjugation 
scheme,  degenerate  four-wave  mixing  in  a saturable  absorber,  was  con- 
ducted on  this  program.  With  company  funding,  we  had  developed  a 
coupled-mode  theory  for  degenerate  four-wave  mixing  in  a two-level 
absorbing  (or  amplifying)  system  of  motionless  atoms.  Here  the  non- 
linear medium  is  a saturable  absorber  with  a resonant  enhancement  at 


48 


each  single  photon  step.  The  solutions  indicate  that  amplified 
reflection,  oscillation,  and  phase  conjugation  occur  even  in  the 
presence  of  strong  absorption.  On  the  present  program,  this  develop- 
ment was  extended  to  include  atomic  motion. 

The  basic  geometry  for  four-wave  mixing  is  given  in  Figure  22. 

The  interaction  consists  of  two  intense  pump  beams  (Ep^,  Ep0)  and  two 
weak  signal  beams  (Eg^,  Eg,,).  The  particular  case  for  which  the  theory 
has  been  developed  considers  all  beams  at  the  same  frequency  and  the 
pumps  and  signals  in  a counterpropagating  geometry  as  shown.  The 
use  of  nondegenerate  frequencies  and  noncollinear  pumps  is  discussed 
in  Section  4.  Physically,  the  intense  beam  Ep^  interferes  with  Eg^ 
to  form  regions  of  high  and  low  intensity.  This  Interference  combined 
with  saturation  in  the  medium  spatially  modulates  the  absorption  and 
index  of  refraction.  This  spatial  modulation  leads  to  a dual-grating, 
or  dual-volume,  hologram  picture,  as  shown  in  Figure  23.  The  forward 
pump  (Epp  and  probe  (Eg^)  form  a small  angle  (large  spacing)  grating, 
while  the  backward  pump  (Ep?)  and  probe  (Eg^)  form  a large  angle  grating. 
In  the  readout  process,  the  backward  pump  (Ep^)  scatters  off  the  small 
angle  grating,  generating  a contribution  to  the  phase  conjugate  signal 
(Eg.,),  and  the  forward  pump  (Ep^)  scatters  off  the  large  angle  grating, 
also  giving  a contribution  to  Eg2-  The  theory  without  atomic  motion 
Indicates  that  the  magnitude  of  the  phase-conjugate  signal  is  a function 
of  three  parameters:  the  line-center  absorption  (gain)  coefficient- 
interaction  length  product,  cx^L;  the  ratio  of  pump  intensity  to  line- 
center  saturation  intensity,  I /Ig^,  w*iere  *So  proportional  to 
(1/T1T2);  and  the  position  in  frequency  of  the  pump  relative  to  the 
line-center  frequency  of  the  saturable  absorber,  y = (to  - u>q)/Am, 
where  Aus  - collision-broadened  linewisth.  For  example,  if  oi  L = 7, 
y = 6,  and  1 ^ / 1 Sq  “ 25,  the  conjugate  return  wave  intensity  is  20X  the 
input  intensity.  if  the  atoms  move  a grating  period  (which  is  of 
order  X)  in  a population  relaxation  time  T^,  the  interference  pattern 
will  he  washed  out,  resulting  in  a reduction  of  the  amplitude  of  the 
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Figure  22.  Four-wave  mixing  scheme 
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conjugate  wave,  Tills  diminution  of  the  conjugate  amplitude  occurs 
even  in  the  case  of  a homogeneously  broadened  lineshape  (1/'1'7  > Av^) . 
This  effect  also  causes  a dependence  of  the  return  wave  amplitude  on 
the  angle  between  the  signal  and  pump  beams.  We  have  developed  (see 
Appendix  C)  and  applied  the  formalism  of  the  "quasi-classical"  density 
matrix  to  compute  the  relevant  nonlinear  couplings,  thus  extending  the 
previous  analysis  to  all  ranges  of  the  atomic  motion  parameter 
(T  = 1/AVp).  The  equations  have  been  solved  exactly  in  third-order 
perturbation  theory  and  in  the  "homogeneous  saturation"  approximation 
(implicit  in  the  stationary  atom  results)  for  arbitrarily  strong  pump 
fields. 

The  solution  for  the  field  amplitude  of  the  phase  conjugate  return 
beam  is  shown  in  Figure  24.  The  figure  gives  the  expected  reduction  of 
the  field  amplitude  compared  to  the  solutions  obtained  previously, 
which  did  not  account  for  atomic  motion.  The  solution  is  a function 
of  T^/T^  and  of  the  angle  between  the  pump  and  signal  beams.  For 
T /T  >>  1 (i.e.,  no  atomic  motion),  the  solution  shows  no  reduction 
over  the  previous  case,  as  expected.  For  angles  greater  than  'v  30° , 
there  is  a general  reduction  in  amplitude  as  T^/T^  decreases  (i.e., 
as  atoms  have  time  to  move  and  wash  out  the  grating).  For  small 
angles,  the  period  of  the  grating  becomes  increasingly  large,  increas- 
ing the  time  it  takes  for  the  atoms  to  move  one  grating  period  and  thus 
minimizing  the  Doppler  motion  reduction.  This  is  exhibited  in  the 
solutions  obtained.  For  small  angles  (such  as  T^/T^  1)»  the  "duc- 

tion  in  amplitude  approaches  a factor  of  two.  This  is  because  the 
large  angle  grating  (smaller  period)  is  becoming  completely  washed  out 
and  can  no  longer  contribute  to  the  return. 

We  have  also  studied  polarization  effects  in  four-wave  mixing  sys- 
tems. The  various  nonlinear  optical  media  in  which  four-wave  mixing 
is  observed  will  exhibit  different  relationships  between  the  polari- 
zations of  the  signal  and  pump  waves  and  the  polarization  and  ampli- 
tude of  the  reflected  (phase  conjugate)  wave.  This  is  because,  even 
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Figure  24 . Reduction  in  amplitude. 
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in  isotropic  materials  such  as  Rases  and  liquids,  the  third-order 
polarizability  contains  two  different  types  o!  tensor  structures,  giving 
rise  to  induced  polarizations  ol  the  foim 


P = 3.  E,  (K 
I s 1 


(K 


(In  crystal  media,  even  more  tensor  structures  can  appear.)  For  example, 
only  the  term  dominates  In  the  single  photon  resonant  process.  In 
the  two-photon  resonant  case,  Bj  probably  dominates.  In  a nonresonantly 
enhanced  process  (CS^,  for  cj  pie),  3^  and  32  can  he  of  comparable 
magnitude . 

Table  1 gives  examples  of  the  expected  output  polarizations, 
relative  efficiency  of  conjugation,  and  angular  dependence  of  the  con- 
jugate signal.  The  examples  are  for  nonresonant,  single-photon  resonant, 
and  two-photoi  resonant  Interactions.  The  arrows  and  dots  indicate  all 
the  nonredundant  directions  ol  the  linear  polarization  states  of  the 
appropriate  beams.  The  efficiencies  Indicate  which  terms  contribute  to 
the  output,  and  the  subscripts  S and  T indicate  single-  and  two-photon 
cases.  For  the  nonresonant  case,  the  first  column  gives  all  beams 
copoiarized.  In  tills  case,  since  the  dot  products  are  all  nonzero,  each 
term  contributes  to  t lie  conjugate  output,  giving  an  efficiency  of 
3 1 + otBj*  hi  tie  second  column,  the  pump  and  signal  are  cross  polarized, 
giving  a zero  dot  product,  but  the  pumps  are  copolarized,  giving  a nonzero 
dot  product  in  the  first  term.  Thus,  the  3j  term  is  the  only  term  remain- 
ing. The  remaining  columns  are  determined  in  the  same  manner.  (The 
absolute  efficiency  is  larger  in  the  resonant  cases.)  The  angular  depend- 
ence arises  in  the  third  columns.  In  this  case,  under  the  assumption  of 
molecular  dipole  scattering,  wacn  the  beams  are  at  90°,  th_  playback  pump 
polarization  Is  oriented  parallel  to  the  probe  propagation  vector,  end 
hence  the  scattered  fields  vanish  In  the  direction  of  the  conjugate 
(.etro  direction).  Thus  operation  is  not  possible  at  this  angle. 
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Table  1.  Effect  of  the  State  of  Polarization  on  the  Conjugate  Output 


perpen  1 ' cular  to  plane. 
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MATERIALS  FOR  FOUR-WAVE  MIXING 


In  this  part,  of  the  program,  we  have  investigated  candidate 
materials  for  four-wave  mixing.  Materials  suited  to  the  blue-green  were 
of  particular  interest,  and  some  are  listed  below. 

The  general  classes  of  nonlinear  materials  useful  for  four-wave 
mixing  can  be  described  under  two  major  categories: 

• Those  that  are  nonresonant  (i.e.,  show  a very  weak 
spectral  dependence  in  the  nonlinearity  x^^that 
results  in  the  phase-conjugate  signal). 

• Those  that  are  resonant  or  show  a pronounced  increase 
in  the  nonlinearity  versus  wavelength,  ns  might  occur 
in  the  vicinity  of  an  allowed  transition  in  a two- 
level  system. 

In  both  cases,  four-wave  mixing  might  occur  as  a result  of: 

• A temporal  modulation  of  the  polarization  at  2ui  that 
scatters  the  signal  wave. 

• A spatial  modulation  of  the  real  part  of  the  refractive 
index  (i.e.,  via  a light- induced  pure  phase  volume, 
hologram) . 

• A spatial  modulation  of  the  absorption  or  gain 
coefficient  (imaginary  part  of  the  refractive  index), 
i.e.,  via  a light-induced  amplitude  hologram. 

• A combination  of  the  above. 

In  any  of  the  three  cases,  the  resonant  media  generally  require  lower 
power  densities  (at  appropriately  chosen  wavelengths)  than  the  non- 
resonant media  for  comparable  reflectivity  of  the  conjugate  wave. 
Examples  of  nonresonant  media  include: 

• Materials  v ith  high  nonlinear  refractive  indices, 
such  as  CS2  (in  the  visible),  germanium  (in  the 
infrared),  and  conjugated  polymers  (such  as 
diacetylenes) . 
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• Materials  that  are  near  a phase  transition  that  may  be 
Induced  local Ly  in  space  by  the  interference  pattern  of 
the  light  field,  such  as  ferroelectr ics  (e.g.,  K1)P) 
near  their  critical  magnetic  fields. 

• Materials  in  which  free-carrier  refractive-index  pro- 
files can  be  easily  Induced,  such  as  semiconductors 
(Si,  CdS,  GaAs,  etc.)  near  and  above  their  band  edges. 

Resonant  media  — both  single-  and  two-photon  transition  media 
might  exhibit  enhancements  of  their  nonlinearity  over  broad  or  ex- 
tremely narrow  spectral  ranges,  depending  on  the  width  of  the  relevant 
resonances  For  laser  pulses  longer  than  a few  nanoseconds,  the  a ingle- 
photon  transitions  appear  more  promising  and  generally  require  extremely 
low  power  densities  for  phase-conjugate  returns  of  a few  percent.  The 
volume  holograms  in  such  media  occur  via  the  nonlinear  saturable  absorp- 
tion of  the  transition,  which  results  In  high  conjugation  efficiencies 
at  power  densities  comparable  to  the  saturation  Intensity  I , Since 
1 Is  inversely  proportional  to  a (the  absorption  cross  section)  and 

Silt 

r (the  recovery  time  for  the  absorption),  phase  conjugation  at  extremely 

2 

low  power  densities  ('t  mW/cm  ) may  be  possible  with  a choice  of  slowly 

recovering  highly  absorbing  media.  An  example  of  such  a material  is 

Na-f luoresceln  (in  orthoborlc  acid),  which  exhibits  remarkable  saturation 

2 

properties  at  optical  power  densities  oi  less  than  500  mW/cm  at  488  tun 

+ 2 

(Ar  laser)  and  at  slightly  higher  powers  ('^  2 W/crn*')  at  501  nm  (HgBr 

laser).  Efficient  (a  few  percent)  backward-wave  generation  has  already 

2 

been  demonstrated  at  relatively  low  power  levels  (“V  1 to  10  kW/cm")  of 
Ar+  laser  radiation  (at  all  lines  from  488  nm  to  5L4.5  nm)  using  a 
1.5-c.m-long  ruby  crystal,  which  acts  as  a broad-band,  three-level, 
blue-green  saturable  absorber  (similar  to  Na-fluorescein  in  orthoborlc 
acid).  Dye  molecules  in  the  fluorescent  state,  as  used  in  liquid  solu- 
tions for  lasers,  also  act  as  excellent  fast-recovery-time  (10  nsec, 
saturab-le  absorbers,  as  might  be  required  for  several  high-dnta-vate, 
large-bandwidth  phase-conjugation  applications.  However,  the  price  for 

speed  is  a higher  operating  power  density  resulting  from  the  correspond- 

2 2 

ingly  larger  saturation  intensities  ('■  100  MW/cm"  to  10  GW/cm"). 


57 


nonlinear  material*  at  2. 7^™'” Ol^laslr ) '’r'ill"1",‘ry  8earch  tor 
lonna,  U10jlM  cs,  „te  promlalng 

appropriate  HP  medium,  absorbing  or  amplifying  could  be  ^ an 

rononant  four-wave  mlxlllg  serial.  ' *”  USUd  38  8 

Certain  specific  examples  of  lmer« 

visible  and  at  2 7 Lim  * 311  candit,ate  materials  in  the 

at  2.7  Mm  are  summarized  in  Table  2a, 

in  the  figure  of  merit  is  tht>  a u A"  lmportant  factor 

is  Particularly  true  for  the  7 eSh°1J  °f  -*  — ial.  This 
efficiencies  are  desired  fromTL  n,Whl.Ch  COnJU8ate  returns  at  high 
the  efficiency  is  generally  n nr<-sonant  material  since  in  this  case 
Power  density.  if  hl  h ffl  °n°tonltdlly  increasing  function  of  the 

much  below  the  damage  threshold  ( ' ^ P°war  densities 

- — nonlinear  ^ 

demonstrated  by  magnifying  optics  that  wi7r  WlU  be 

SUch  beams  (to  optimal  levels)  at  th  ^ ^ ^ denslty  of 

■revels)  at  tn«  phase  conjugates 

n*  THREE-WAVE  MIXING 

tn  addition  to  SBS  and  four-wave  mixing 
using  a parametric  difference  fr  3 ^ conjugatlon  scheme 

examined.  This  process^  £1  * ~l0n  P—  - briefly 

-2—  scheme  "L“  ^ ^ 

multimode  flbe  ptlc  waveguide  used  f„r  " dl',P"Si°"  lntroduc«1  in  a 

CeSS  Mk'S  “8=  '■»  quadratic  ^ Pr°' 

ceptlbint,  present  In  8"8~ 

Figu«  ’I(“‘tlI„0Ltl,tl,lTC  "h<,ae'CO"JUSatlc"'  Is  shown  ln 

18  •*»'“  « nonlinear  0^™'.^"'’,''  ““Ve  f Ce')"‘’”C>'  u3 

-8  .-ing.  a new  w.  (UuT>  ""  “ V 

“2(“2  ■ “3  - “,)  I*  Senerated  that  is  the  ere"Ce  fre''“e"'>' 

»uve.  That  is,  the  „,lxl  P 88  “njugate  of  the  signal 

■■  - 


58 


7182  14 


lw)l 

(<) 


(b) 


Figure  25.  Three-wave  mixing  schemes. 


Table  2.  Candidate  Materials 


Laser  (Wavelength) 

Material 

Typical  Power 
Densities'1 

Dye  lasers  (e.g., 

Dye  solutions  (e.g., 

V100  MW  to  1 GW/cm2 

Rhodamine  6G  at 

10“^M  Oxazine  in 

590  nm) 

ethanol) 

Doubled-YAG  (532  nm) 

Ruby  (0.05%) 

^1  to  10  kW/cm2 

Si  (p  > 100  (2-cm) 

2 

Vl  MW/ cm 

HgBr  (501  nm) 

Na- fluorescein 

SI  W/cm2 

Ruby  (0.05%) 

^1  to  10  kW/cm2 

Ar+  (488  nm/ 

514.5  nm) 

Na-f luorescein 

Si  W/cm2 

Ruby  (0.05%) 

'vl  to  10  kW/cm 

Ruby  (694.3  nm) 

Rubrocyanine  in 
methanol 

Vl  kW/cnT 

CdSe , CdS 

VI  to  10  MW/cm" 

Cryptocyanine  in 
methanol 

^•50  MW/cm2 

For  conjugate  returns 

of  a few  percent. 

other.  Thus,  if  E^  is  the  input  signal  wave,  a wave  emerges  at  the 
exit  of  the  medium  where  E^  = KE*.  If  the  pump  wave  E^  is  an  intense 
field,  the  wave  E^  may  exhibit  amplification  (K  > 1).  This  is  a result 
of  conservation  of  energy.  Here  the  energy  of  the  pump  wave  is  coupled 
through  the  polarization  generated  at  frequency  ^ by  the  non- 

linear mixing  to  the  field  E^.  If  the  pump  wave  frequency  is  set 
equal  to  2oo^,  then  the  down-converted  wave,  will  have  the  same  fre- 
quency as  the.  input  wave  and  will  be  the  phase  conjugate.  The  latter 


60 


process  (shown  in  Figure  25(b))  is  of  primary  interest  for  most  applications. 
From  an  examination  of  the  solution  of  the  differential  equations  given  by 
Yariv  describing  three-wave  mixing,  one  observes  that  for  multimode 
fields  (general  aberrated  beam),  each  mode  of  the  signal  field,  E^(w^), 
interacts  only  with  its  counterpart  mode  at  o^.  This  is  a result  of 
the  need  to  "phase  match"  in  the  transverse  direction  (across  the  beam) 
and  in  the  propagation  direction.  This  "three-dimensional"  phase  match- 
ing will  ensure  exact  phase  conjugation  . The  need  to  phase  match  each 
mode  to  achieve  exact  phase  conjugation  is  a serious  limitation  on  three- 
wave  mixing.  This  implies  that  only  a narrow  bundle  of  rays  or  a few 
modes  that  are  nearly  phase  matched  on  alignment  of  the  beams  will  have 
sufficient  amplitude  to  be  observed  and  are  the  only  modes  phase 
conjugated. 

Since  the  three  waves  in  this  example  are  traveling  in  the  same 
direction  through  the  crystal,  there  exists  a need  to  be  able  to  dis- 
tinguish the  input  wave  from  the  conjugate  wave.  This  is  done  by  making 
the  polarization  of  the  generated  wave  * KE^*  orthogonal  to  the  signal 
wave  E^  (waves  E^  and  E^  have  the  same  polarization).  This  requires 
a material  that  satisfies  the  type  II  phase-matching  condition.  Such 
crystals  might  be  Te,  Ag^AsS^,  CdSe,  or  KDP. 

An  experimental  demonstration  of  phase  conjugation  by  three-wave 
mixing  was  done  by  Avizonis  et  al.  (Appl.  Phys.  Lett.  Vol.  31,  pp.  435- 
437,  Oct.  1977)  at  1.06  pm.  The  crystal  used  was  lithium  formate. 

Conjugate  signals  0.5%  of  the  input  signal  were  observed  to  be  limited 
by  the  ability  to  phase  match  (Ak  = k^  “ k^  - 1 0) . These  experiments 

showed  that  the  conjugate  signal  observed  was  an  extremely  sensitive 
function  of  the  alignment. 

Figure  26  shows  a conceptual  realization  of  this  process  for  a 
HEL  using  a shared  aperture.  An  illumination  pulse  is  initiated  in  the 
same  manner  as  in  Figure  1 to  obtain  a glint  return.  The  glint  return 
passes  through  the  parametric  amplifier  to  generate  the  conjugate  wave- 
form, which  is  selected  by  the  polarizer.  This  wave  is  further  amplified 
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Figure  26.  Phase  conjugation  scenario  using 
three-wave  mixing. 


and  sent  back  through  the  power  amplifier  to  the  target.  One  advantage 
of  this  process  is  the  free  choice  of  the  frequency  shift  of  the  con- 
jugate wave  with  respect  to  the  input  wave.  A major  disadvantage  of 
this  scheme  is  that  the  path  containing  the  Faraday  rotator  isolator 
is  not  included  in  the  conjugation  process  and  is  not  corrected. 
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SECTION  4 


SYSTEMS  APPLICATIONS 

A.  SYSTEMS  OVERVIEW  AND  SCOPE 

We  have  tentatively  examined  several  possible  applications  for 
nonlinear  phase  conjugation  (NPC)  systems.  These  applications  fall 
basically  into  three  categories:  (1)  local  corrections  for  optical 
defiguring  or  local  aberrating  paths  (e.g.,  laser  media),  (2)  relay 
mirrors  for  path  distortions  up  to  but  not  including  the  relay  mirror, 
and  (3)  target-referencing  systems  (for  all  classes  of  distortion). 

The  basic  objectives  are  (1)  to  identify  any  systems  limitations  that 
may  occur  because  of  the  inherent  nature  of  NPC  and  (2)  to  identify 
new  techniques  for  using  NPC  that  could  have  a major  Impact  on  systems 
design. 

We  have  emphasized  operation  at  visible  wavelengths  for  two 
reasons.  First,  assuming  that  the  required  laser  development  is  forth- 
coming, visible  operation  offers  a major  potential  for  improved  overall 
system  operation  In  future  long-range  systems.  More  explicitly,  ex- 
cepting turbulence-induced  problems,  the  shorter  wavelengths  both  offer 
reduced  diffraction  spreading  and  appear  to  be  less  susceptible  to 
atmospherically  Induced  nonlinear  effects  because  of  lower  absorption. 
Second,  it  is  only  at  visible  wavelengths  that  the  path  compensation 
advantages  of  NPC  systems  strongly  dominate  those  of  conventional 
adaptive-optic  systems.  Specifically,  a 4-m  aperture  operating  at 
visible  wavelengths  through  the  atmosphere  will  require  about  3,000 
control  systems  and  deformable  mirror  actuators  or  control  points. 

Hence,  NPC  offers  the  potential  of  replacing  a very  complex  system 
with  a single  NPC  correcting  device. 

In  summary,  we  believe  that  operation  at  visible  wavelengths  offers 
major  advantages  for  ground-to-space  missions,  and  that  a large  NPC 
system  offers  potential  for  compensating  the  substantial  loss  in  per- 
formance produced  by  atmospheric  turbulence  and  its  secondary  effects 
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preceding  page 


(amplitude  scintillation)  that  occur  at  visible  wavelengths.  That  is, 
we  believe  that  amplitude  preservation  (of  the  reference  wave  field)  in 
the  retransmitted  wavefront  may  be  essential  for  optimium  correction  at 
visible  wavelengths  and  that  NPC  offers  a promising  way  of  achieving  it. 

Although  we  have  high  confidence,  based  on  our  experimental  work, 
that  NPC  will  achieve  excellent  compensation  for  laser  beams  that  share 
a fully  common  path  with  the  reference  (or  probe)  wave,  such  a path  is 
not  easy  to  achieve  in  operational  systems.  For  example,  the  point 
ahead  angles  required  when  operating  against  orbiting  targets  (although 
only  tens  of  microradians)  appear  to  cause  sufficient  lack  of  path 
commonality  to  degrade  system  performance  badly.  This  effect,  commonly 
called  the  point-ahead  isop.lanatic  problem,  is  more  fully  discussed  in 
Appendix  D.  Since  no  solution  to  this  problem  (at  visible  wavelengths) 
appears  to  exist,  we  have  excluded  this  class  of  system  from  explicit 
consideration  in  this  report. 

The  relay  systems  offer  appreciable  advantages  in  this  respect. 
Specifically,  one  can  position  or  fly  a reference  ahead  of  the  relay* 
and  thereby,  in  principle,  fire  the  laser  back  over  substantially  the 
identical  propagation  path  taken  by  the  reference.**  However,  practical 
considerations  associated  with  reference  size  motivate  consideration  of 
operational  modes  that  do  not  employ  a fully  common  laser  and  reference 
path,  as  discussed  below. 

Although  isoplanatic  effects  provide  the  key  argument  against  long- 
range  target  referencing  systems,  there  are  others.  In  many  respects, 
the  NPC  system  shares  the  same  problems  as  the  conventional  phase- 
conjugate  adaptive  optical  systems:  (1)  convergence  time  is  limited  by 


This  concept  (as  far  as  is  known)  was  independently  conceived  and  pro- 
posed to  NASA  by  Hughes  and  Rockwell  International  in  1974  proposals  to 
transmit  laser  power  from  a ground  station  to  a cooperative  satellite. 

The  surface  and  associated  air  mass  of  the  Earth  move  a bit  during  the 
transit  time  from  the  top  of  the  atmosphere  to  the  surface  and  back, 
but  this  produces  only  a minor  perturbation  from  path  commonality. 
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(typically*)  four  or  five  round-trip  propagation  times  and  (2)  with 
static  targets,  the  (nearly  diffraction  limited)  focused  beam  forms  on 
the  dominant  glint  or  highlight  as  a reference,  a reference  that  may 
"burn  off"  as  operation  proceeds,  shifting  the  aim  point  to  the  next 
glint,  etc.  The  generally  proposed  solution  to  this  latter  problem  is 
to  offset  the  laser  beam  from  the  reference  by  electronic  override 
techniques  (such  techniques  have  been  experimentally  demonstrated).  At 
the  beginning  of  this  contract,  it  was  not  self-evident  whether  such  an 
offset  could  be  achieved  in  an  NPC  system  since  its  natural  operation 
is  to  focus  “he  laser  beam  directly  on  the  dominant  glint,  and  simple 
electronic  tilt  override  is  not  possible.  We  have  found,  however,  that 
it  is  possible  to  achieve  an  electrooptic  tilt  override  by  several 
techniques.  One  approach  is  to  employ  a four-wave  mixer  NPC  by  perturb- 
ing the  pump  angles  and/or  frequencies  somewhat  (this  is  discussed  below). 

Another  problem  that  carries  over  from  conventional  phase  conjuga- 
tion systems  is  target-motion-inuuced  Doppler  shifts.  These  are  gener- 
ally removed  in  a conventional  adaptive  optical  system  by  either  a 
Doppler  tracking  local  oscillator  or  electronically  by  a second  (track- 
ing) local  oscillator  downconverting  the  i.f.  output  of  the  optical 
heterodyne  detector.  Again,  we  have  found  that  appropriate  frequency 
and  angle  offsets  of  the  pumps  can  achieve  a frequency  shift  adequate 
for  Doppler  compensation.  In  fact,  it  is  possible  to  achieve  simultan- 
eous pointing  offset  and  Doppler  compensation  and  yet  maintain  perfect 
phase  matching  in  the  four-wave  mixer.  (This  is  also  discussed  below.) 
Thus,  although  we  do  not  explicitly  consider  target  referencing  systems 
in  this  report,  we  do  consider  two  of  the  key  problems  associated  with 
them  since  they  share  these  problems  with  the  relay  systems. 

System  applications  to  local  correction  systems  can  be  further 
broken  down  into  the  following  categories  for  more  detailed  examination: 


Assuming  that  multiwavelength  conjugation  is  feasible,  multi wavelength 
operation  decreases  convergence  time.  However,  if  the  amplitude  dis- 
tribution in  the  reference  return  wavefront  must  be  preserved,  this 
will  considerably  increase  convergence  time. 
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(la)  oscillator  compensation,  (lb)  MOPA  compensation,  (lc)  satellite- 
based  optical  train  compensation.  Likewise,  several  approaches  to 
relay  systems  applications  are  considered.  These  systems  employ  NPC 
compensation  for  a ground-based  system  operating  to  an  orbiting  relay. 
The  scope  of  the  attempted  wavefront  correction  progressively  increases 
through  this  sequence.  In  applications  la  and  lb,  the  aim  is  primarily 
to  correct  for  the  laser  medium  distortions  and,  in  some  cases,  for 
some  of  the  intracavity  optical  elements.  In  application  lc,  the  key 
compensation  is  for  an  orbiting  primary  mirror  and  the  lasing  medium. 

In  relay  applications,  the  aim  is  to  compensate  the  entire  ground-based 
optical  train,  including  the  laser  medium  and  the  atmospheric  portion 
of  the  propagating  path. 

We  have  included  several  structurally  different  approaches  in 
each  application;  in  particular,  we  have  attempted  to  examine  both 
oscillator  and  MOPA  alternatives  when  they  could  be  found.  Further, 
we  have  included  systems  appropriate  for  both  short  pulse  and  cw  (or 
long  pulse)  operation  to  best  accommodate  future  development  in  laser 
technology.  Although  our  major  emphasis  is  on  visible  systems,  we 
believe  that  the  basic  concepts  are  applicable  at  any  wavelength. 

While  investigating  the  relay  systems,  we  found  that  it  is  ex- 
tremely desirable  (perhaps  essential)  to  provide  both  a tilt  override 
capability  (as  previously  discussed  for  target  referencing)  and  a focus 
override  (to  prevent  the  conjugated  wavefront  from  being  simply  retro- 
reflected  with  a convergence  equal  to  the  divergence  of  the  incident 
beam).  This  requirement  comes  about  because  the  geometry  of  the  trans- 
mitted beam  and  the  reference  beam  may  not  be  the  same.  In  particular, 
it  is  decidedly  convenient  to  use  a reference  that  is  much  smaller  than 
the  relay  mirror.  At  the  extreme,  we  could  have  a situation  in  which 
the  transmitted  beam  was  essentially  collimated  and  the  reference  beam 
was  generated  by  a point  source  located  (by  approximately  the  point- 
ahead  angle)  ahead  of  the  relay  mirror.  The  geometrical  difference  in 
the  paths  taken  by  these  two  waves  leads  to  both  a focus  offset  require- 
ment in  the  NPC  system  and  to  additional  isoplanatic  problems  not  pres- 
ent when  the  transmitted  and  reference  waves  have  similar  geometries. 
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B. 


SYSTEMS  INVESTIGATED 


1.  Local  Correction  Systems 

A limited,  but  nevertheless  very  important,  application  of  non- 
linear phase  conjugation  is  to  the  internal  compensation  of  oscillators 
and  master-oscillator  power-amplifier  (MOPA)  systems.  The  potential  ad- 
vantages over  conventional  intracavity  adaptive-optics  systems  are  (1) 
simultaneous  correction  on  N laser  lines,  each  with  its  own  distortion, 
appears  feasible;  (2)  much  higher  spatial-frequency  distortions  can  he 
corrected;  (3)  very  much  higher  bandwidth  corrections  can  be  a 'loved; 
and  f4)  ultimate  cost  should  be  lower  because  the  system  is  intrinsi- 
cally simple. 

a . Oscillator  Compensation 

Opj  basic  approach  to  oscillator  compensation  is  illustrated 
in  Figure  27.  rn  a general  sense,  three  key  elements  are  included: 

(1)  a plane  wavefront  generator,  (2)  a lasing  medium,  and  (3)  a non- 
linear phase  conjugator.  In  the  system  illustrated,  the  plane  wavefront 
is  produced  by  a pinhole  spatial  filter,  and  the  phase  conjugator  is  an 
SBS  system.  To  facilitate  the  discussion  of  these  systems,  * will  use 
transmission  optics  and  low-power  implementations.  In  some  cases,  we 
will  also  discuss  implementations  and  problems  appropriate  to  high-power 
applications . 

There  are  several  fairly  basic  problems  with  the  system  depicted  in 
Figure  2/.  First,  there  is  a problem  with  pinhole  burnout.  Second, 
because  the  pinhole  transmission  with  badly  distorted  wavefronts  may, 
during  Initial  start-up,  reduce  the  ensemble  gain  below  SBS  threshold 
conditions,  th'  re  is  a problem  with  threshold  start-up.  Third,  there 
are.  two  problems  with  the  beam  splitter:  a large  portion  of  the  plane 
wavefront  returning  from  the  pinhole  is  lost*  by  the  splitter,  and, 

* 

The  lost  energy  can  of  course  be  largely  recovered  by  adding  a mirror 
which  redirects  the  beam  back  in  the  direction  of  the  output  wavefront. 
Such  a recovery  system  must  act  like  an  interferometer  so  the  correct 
phasing  and  angular  orientations  must  be  carefully  maintained. 
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Figure  27.  A basic  approach  to  oscillator  compensation.  The 
beam  splitter  would  be  replaced  by  a diffraction 
grating  at  high  power  levels  and  the  illustrated 
splitter  ratio  gives  about  1035  feedback. 
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Figure  28.  A more  efficient  approach  to  oscillator  compensa- 
tion. By  employing  a phase  conjugator  which  rotates 
the  plane  of  polarization  and  a polarization- 
sensitive  beam  splitter,  the  power  lost  by  the 
splitter  is  largely  avoided  and  ^he  power  density 
on  the  pinhole  is  reduced  by  a factor  of  three 
(for  the  same  feedback). 
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consequently,  the  power  density  on  the  pinhole  is  higher  than  need  be. 
Fourth,  the  frequency  offset  produced  in  a basic  SBS  conjugator  pro- 
gressively accumulates  and  eventually  walks  off  the  gain  profile  line 
width  of  the  lasing  medium. 

In  Figure  28,  we  have  corrected  some  of  these  defects  by  incor- 
porating a polarization  switching  phase  conjugator  using  the  polarization 
rotation  schemes  discussed  in  Section  3.  Thus,  the  wavefronts  reflected 
off  the  conjugator  are  polarized  parallel  to  the  splitter  (normal  to 
the  plane  of  the  paper),  and  the  polarization  of  the  wavefront  trans- 
mitted back  via  the  pinhole  system  is  rotated  by  90°  by  a double  pass 
through  a quarter-wave  plate.  The  polarization  plate  is  chosen  to  pass 
this  polarization  substantially  without  reflection.  For  example,  it  may 
be  an  interference  filter  operated  near  Brewster's  angle  for  the  filter, 
or  it  may  be  a diffraction  grating  chosen  to  give  high  efficiency  at  this 
polarization.  Thus,  for  10%  feedback,  for  example,  the  beam  splitter 
can  be  designed  to  transmit  10%  of  its  incident  energy  to  the  pinhole, 
as  contrasted  to  the  33%  transmitted  by  the  system  shown  in  Figure  27. 
Further,  the  "lost"  energy  from  the  second  pass  of  the  splitter  in 
Figure  27  is  eliminated.  To  avoid  the  frequency  walk  off  problem,  we 
have  replaced  the  SBS  conjugator  with  a four-wave  conjugator  in 
Figure  28. 

Figure  29  shows  a pump  supply  to  the  system  shown  in  Figure  27 
wherein  the  sustained  pump  is  achieved  via  oscillator  feedback.  Since 
the  pump  supply  will  not  function  during  initial  start-up,  we  have  added 
a one-pulse  pulsed  oscillator,  which  supplies  the  required  pump  power 
during  the  start-up.  We  have  taken  advantage  of  the  degeneracy  of  the 
four-wave  mixing  system  to  isolate  the  two  types  of  pump  by  angular 
separat ion. 

In  addition  to  the  technology  problems  presented  by  the  phase  con- 
jugator, the  pinhole,  and  the  splitter  elements,  there  are  basic  problems 
associated  with  the  confocal  resonator  configuration  associated  with 
made  volume  filling. 


71 


**67  -7 


Figure  30.  The  "rotary" 
Pinhole  system.  A pair 

an  pr°ta.tln8  wheels  with 
an  annuiar  hole  in  each 
wneel  forms  a pinhole 
with  rotary  edges.  in 

Dinh  r*  heatln8  at  the 
Pinhole  edges  ls  distri- 
buted over  an  area 
hundreds  of  times  reater 

pinhole 3 
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Figure  31.  A "figure  eight"  ring  laser  oscillator  with  a phase 
conjugator  for  gain  medium  compensation.  The  gain 
medium  is  double  passed  in  the  same  direction  in  two 
differing  polarization  states  before  exiting  the 
output  polarization  separation  systems. 
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We  have  explored  alternatives  (such  as  Fabry-P^rot  resonators)  to 
the  pinhole  spatial  filters  without  finding  anything  of  strong  promise. 
Currently,  it  appears  that  t'->e  best  approach  would  be  to  use  rotary 
metal  disks  to  distribute  the  pinhole  peripheral  heating  over  a much 
wider  area,  as  illustrated  in  Figure  30.  The  pinhole  system  should  be 
somewhat  canted  such  that  reflected  energy  from  the  periphery  of  the 
pinhole  does  not  re-enter  the  oscillator. 

We  have  also  explored  several  possible  ring-resonator  oscillator 
configurations;  the  most  promising  one  to  date  is  illustrated  in  Fig- 
ure 31.  Basically,  the  feedback  follows  a figure  eight  path.  The  upper 
loop  of  the  "8"  is  the  usual  optical  feedback  loop,  while  the  lower  loop 
is  a double-pass  loop  in  which  the  polarization  state  of  the  wavefront 
transmitting  the  amplifier  is  first  normal  to  the  plane  of  the  loop. 

The  wavefront  reflects  off  the  output  polarization  filter  to  the  phase 
conjugator,  reflects  as  a conjugated  wavefront,  has  its  polarization 
state  rotated  to  the  plane  of  the  loop,  passes  the  amplifier  once  more, 
and  then  exits  through  the  exit  polarizer.  Two  optically  active  rotators 
are  required  to  set  up  the  double  pass  of  the  lasing  media.  Note  that 
the  phase  conjugator  is  operating  off-axis  in  a mode  contrary  to  its 
natural  functioning.  This  operation  can  be  obtained  in  a four-wave 
conjugator,  as  described  in  Appendix  E, 

Another  approach  to  the  "plane-wave-make"  element  of  the  basic 
oscillator  system  is  illustrated  in  Figure  32.  This  approach  employs 
a four-wave  mixer  with  pumps  derived  from  the  probe  or  signal  wavefront 
as  was  the  case  in  the  oscillator  system  of  Figure  29.  In  the  present 
case,  however,  the  distortions  present  in  the  probe  wavefront  are  retained 
on  the  pump.  For  example,  consider  a wavefront  distortion  that  is  essen- 
tially a one-dimensional  distortion  with  variation  in  the  y direction. 

This  distortion  maps  into  a z-direction  distortion  in  the  pump  field, 
thereby  destroying  the  regular  spacing  in  the  gratings  that  reflect  the 
probe  wave.  Thus,  the  grating  efficiency  is  reduced  by  the  dephasing  of 
the  scattered  probe  signal,  and  the  reflected  signal  is  accordingly 
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reduced.  Accordingly,  the  single-pass  gain  for  undistorted  wavefronts 
is  much  greater  than  for  distorted  wavefronts,  and  the  system  will 
oscillate  with  modes  that  correspond  to  such  wavefronts.  Distortions  in 
the  orthogonal  direction  can  produce  corresponding  reductions  in  gain  by 
operating  in  tandem  off  a second  four-wave  mixer  whose  pump  direction 
is  normal  to  the  original  set. 

b . MOPA  Compensation 

Because  of  the  double-pass  nature  of  phase  conjugation  correc- 
tion systems,  we  have  found  that  amplifiers  are  more  natural  candidates 
for  correction  than  are  oscillators.  One  approach  to  MOPA  compensation 
is  illustrated  in  Figure  33.  In  this  figure  the  oscillator  and  amplifier 
are  separately  compensated.  The  oscillator  system  is  basically  the 
same  as  previously  illustrated  in  Figures  27  and  28  with  the  mirror 
and  the  four-wave  mixing  conjugator  forming  the  resonator  system.  A 
portion  of  the  oscillator's  circulating  power  is  split  out  by  splitter 
Sj,  passes  through  the  amplifier,  is  conjugated  and  rotated  in  polariza- 
tion on  reflection,  and  again  passes  through  the  amplifier,  where  the 
power  level  is  increased  and  the  amplifier  distortions  cancel  in  the 
second  (output)  pass  through  the  amplifier.  In  order  to  improve 
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Figure  32.  Oscillator  with  a nonlinear  plane-wave 
generator.  Functionally  the  four-wave 
mixer  substitutes  for  the  pinhole  element 
of  Figure  27. 
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Figure  33.  A MOPA  system  with  separate  phase  conjugators  to  control  the 
wavefront  of  the  oscillator  and  the  wavefront  of  the  ampli- 
fier output. 
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coupling  efficiency.  The  splitter  is  a polarization  filter  designed 
to  pass  substantially  100%  of  the  wave  that  impinges  on  it  polarized 
in  the  plane  of  the  paper.  Thus,  polarization  isolation  minimizes  the 
likelihood  that  feedback  from  the  amplifier  will  capture  the  oscillator 
system.  Note  that  the  amplifier-compensating  four-wave  mixer  is 
arranged  in  such  a way  that  its  pump  fields  are  supplied  by  the  circu- 
lating fields  within  the  oscillator. 

Many  other  MOPA  configurations  are  feasible.  However,  since  MOPAs 
are  used  exclusively  in  the  systems  that  follow,  we  defer  further  dis- 
cussion to  the  discussion  of  those  systems. 

c . Satellite  Optical  Train  Compensation 

This  general  systems  concept  is  illustrated  in  Figure  34. 
Basically,  this  is  a MOPA  system  in  which  the  oscillator  is  injecLed 
into  the  power  amplifier  via  a low-efficiency  diffraction  grating  coupler 
embossed  on  the  primary  mirror.  In  this  way,  both  the  primary  mirror 
and  secondary  mirror  (and  additional  optical-train  elements  (not  illus- 
trated) as  well  as  the  gain  medium  of  the  amplifier  are  double  passed. 
Thus,  all  these  elements  have  their  OPD  distortions  (largely)  compen- 
sated by  the  familiar  double-pass  operation  of  the  phase  conjugator. 

In  more  detail,  the  oscillator  wavefront,  which  we  assume  to  be  of 
high  quality,  is  focused*  to  a spot,  most  likely  in  the  obscuration  hole 
of  the  secondary,  and  diverges  to  fill  the  primary.  A nominal  percentage, 
say  1%,  of  this  wavefront  is  diffracted  off  a shallow  diffraction  grating 
to  a direction  that  is  collimated  with  that  of  an  undistorted  feed  source 
for  the  telescope.  The  remainder  reflects  in  the  zero  order  of  the 
grating  in  the  output  direction.  The  grating  is  circularly  symmetric, 
and  its  period  is  a function  of  radius  because  the.  required  diffraction 
angles  are  a function  of  radius. 

The  diffracted  component  exits  the  telescope  as  a collimated  beam 
if  the  primary  and  secondary  are  well  figured,  traverses  the  amplifier 


We  have  examined  the  feasibility  of  focusing  to  a diffraction-limited 
spot  in  a space  environment  on  other  programs.  It  appears  that  this 
can  be  accomplished  without  breakdown. 
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Satellite-based  system.  Although  the  power  levels  and  flows  are  illustrated 
from  right  to  left  in  the  lower  half  and  reversed  above  it  will  be  understooi 
that  the  actual  power  flows  are  two  way  over  both  paths. 


and  thereby  pick,  up  an,  media  path  distortions.  Is  conjugate-re  letted, 
again  passes  through  the  amplifier,  and  exits  It  with  the  ampin  er 

distortions  stripped  off.  The  majority  (99%,  of  the  output  power  is 

reflected  in  the  aero  order  off  of  the  primary  mirror  and  exits  the 

system  as  a high-quality  wavefront. 

If  the  primary  is  distorted,  then  the  double-pass  system  w 
result  in  a very  substantial  primary  compensation.  However,  with  very 
purge  distortions,  the  angular  differences  between  the  reflected  and 
injected  wavefronts  cause  some  deterioration  In  performance. 

Although  we  have  greatly  relieved  much  of  the  requirement  for 
precision  optical  figuring  on  the  mnin  telescope,  the  price  has  been  a 
requirement  for  precision  Injection  of  the  oscillator  reference.  More 
specifically,  any  distortions  in  the  injected  wavefront  impinging  on 
the  primary  are  directly  mapped  onto  the  output  wavefront.  For  example, 
if  the  injection  focus  is  shifted  a distance  D.  as  Illustrated  in 
figure  35,  a change  In  optical  path  difference  fro.  the  outer  edge  »t 
the  mirror  to  the  inner  (obscuration)  edge  -111  be  produced;  from  the 
geometry  of  Figure  35.  this  change  can  be  shown  to  be 


AOPD 


Mcos  °min  - cos  °max) 


ave  the  angles  measured  to  these  two  edges,  as 


where  °max  and  ^mih  , ,n 

illustrated  in  the  figure.  For  typical  values  (0 max  - 

this  gives 


6 . =6°), 
min 


AO PH  -10.13)0  . 


Thus  0 must  be  held  to  a tolerance*  of  about  plus  or  minus  one  wavelength 
to  hold  this  class  of  induced  error  within  reasonable  tolerances,  t.  s 
Is  not  an  easy  objective  to  achieve.  Although  this  is  only  a focus 


*This  tolerance  is  not  unique  to  NPC-type  systems, 
problem  with  all  local  referencing  wavefront  erro 

systems. 


it  is  a typical 

-sensing/ compensation 
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error  to  first  order,  it  is  not  clear  how  one  could  obtain  the  information 
to  compensate  it. 

A small  amount  (nominally  i%)  of  the  power  radiated  from  the  pri- 
mary is  diffracted  back  into  the  oscillator  system,  where  (at  higher 
levels)  it  presents  a problem  in  its  potential  for  capturing  the 
oscillator.  It  is  this  feedback  problem  that  has  forced  us  to  employ 
a low-efficiency  grating  and  concomitant  high-power  oscillators.  We 
discuss  more  sophisticated  approaches  to  feedback  suppression  below. 

If  high  oscillator  power  levels  must  be  used,  then  the  nominal 
oscillator  of  Figure  34  may  become  a MOPA  system  (as  in  Figure  33,  for 
example).  One  way  of  implementing  an  intermediate  power  amplifier  is 
to  "reflex"  or  iterate  the  basic  concept,  as  illustrated  in  Figure  36. 

Another  method  of  attacking  the  feedback  problem  is  to  use  a 
quarter-wave  plate  ahead  of  the  conjugator  to  rotate  the  plane  of 
polarization  in  the  second  pass.  The  grating  coupling  could  now  be 
increased  to  5 to  7%  with  a decrease  in  required  driver  power  by  a 
factor  of  five  to  seven.  The  5 to  7%  feedback  power  is  suppressed  by  a 
polarization  filter  between  the  primary  and  the  oscillator. 

2 . Orbiting  Relay  Systems 

a . Basic  System  Configurations 

Basically,  we  have  considered  three  classes  of  satellite  sys- 
tems in  which  we  are  transmitting  power  to  a satellite  relay,  as  illus- 
trated in  Figure  37.  In  all  three  systems,  a four-wave  mixer  called  <J>* 
is  used  to  compensate  the  propagation  path  distortions  introduced  by  the 
atmosphere,  the  ground-based  optical  train,  and  an  aberrating  laser 
medium. 

All  three  systems  employ  MOPA  laser  configurations,  in  which  a well- 
controlled  wavefront  (generated  by  a reference  oscillator)  is  used  to 
probe  the  path  distortions.  For  the  moment,  we  take  the  satellite  to 
be  frozen  in  orbit  for  conceptual  purposes.  System  A is  a cw  or  long- 
pulse  system  that  locates  the  oscillator  on  the  satellite.  The  reference 
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oscillator  output  wavefront,  which  is  assumed  to  be  i.ndistorted, 
transmits  downwards  "through"  the  relay  mirror,  picks  up  (probes)  the 
atmospheric  and  amplifier  distortions,  "reflects"  from  the  phase  conju- 
gator,  is  again  amplified,  retransmits  through  the  beam  director  and 
the  atmosphere,  and  exits  the  atmosphere  as  an  undistorted  beam. 

System  B employes  a ground-based  reference  and  injects  a train  of 
short  pulses  (nominally  about  1 usee)  at  a desired  pulse  repetition 
frequency.  Each  injected  oscillator  pulse  generates  a triplet  of  pulses 
(as  illustrated  in  Figure  37),  spanning  (roughly)  the  round-trip  propa- 
gation time  in  duration.  The  first  pulse  of  the  triplet  reflects  off 
a corner  reflector  or  other  retro  system,  returns  after  the  round  trip 
delay  as  a second  pulse  with  superimposed  atmospheric/amplifier  distor- 
tions, is  conjugated  and  reamplified,  and  returns  to  the  satellite  relay 
mirror  as  the  higher-power  and  wavefront-compensated  third  pulse  of  the 
criplet.  The  conjugator  or  amplifier  is  then  temporarily  deactivated 
such  that  no  furuier  iterations  occur  before  the  next  oscillator  pulse 
is  injected  into  the  system. 

System  C is  similar  except  that  a longer  oscillator  pulse  (equal 
to  the  round  trip  time)  is  employed  and  the  oscillator  is  shut  off 
after  the  initiating  pulse.  The  system  gain  is  controlled  such  that  the 
system  iterates  indefinitely  on  the  initial  pulse,  giving  a system 
operation  that  is  effectively  cw  (i.e.,  a system  that  transmits  contin- 
uously for  several  seconds  or  more).  For  conceptual  purposes,  we  have 
shown  the  retro-reflector  in  this  case  as  a number  of  mini-retros 
embedded  in  the  relay  mirror;  however,  in  fact  it  would  be  a more  con- 
ventional retro-reflector  located  in  front  of  the  relays,  as  discussed 
below. 

We  have  identified  three  basic  classes  of  functional  problem  that 
occur  in  all  of  these  systems.  First,  during  the  round  trip  propagation 
time,  the  reference  moves  from  its  old  position  in  such  a way  that  the 
system  points  behind.  That  is.  It  refocuses  back  at  the  old,  rather 
than  the  current,  reference  retro  or  oscillator  position.  Second,  each 
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received  pulse  is  generally  Doppler  shifted  in  frequency  from  that  of 
the  oscillator,  as  illustrated  in  Figure  38.  In  fact,  with  the  System 
C,  the  Doppler  shift  progressively  accumulates  such  that  the  return 
pulses  will  eventually  have  their  line  center  shifted  off  of  the  gain 
line  profile,  and  the  system  will  cease  to  amplify  (and  to  operate). 
Third,  with  small  references,  the  conjugator  system  refocuses  to  the 
dimensions  of  the  retro  or  reference  oscillator,  or  as  closely  as  it 
can  subject  to  diffraction  limitations,  whereas  the  desired  focus  on  die 
laser  beam  corresponds  to  a full  illumination  of  the  relay  mirror. 

The  best  solution  to  the  "point-behind"  problem  is  to  locate  the 
reference  oscillator  or  retro  system  in  front  of  the  satellite  relay 
mirror  by  a distance  corresponding  to  the  point-ahead  angle  (thereby 
minimizing  the  isoplanatic  problems  associated  with  point-ahead  systems). 
Further,  it  avoids  the  problem  of  having  to  operate  "through"  or 
"around"  the  relay  mirror.  We  will  call  this  the  "carrot-on-the-sL iek" 
approach,  as  illustrated  in  Figure  39.  The  required  angular  separation 
between  the  reference  and  the  relay  (centers)  is 

2v 

bp  * — COS  D , (1) 

where  v^  is  the  orbiting  tangential  velocity,  c is  the  velocity  of  light, 
and  0 is  the  azimuthal  angle.  Unfortunately,  since  the  required  exact 
spacing  is  a function  of  the  azimuthal  angle,  it  varies  with  overhead 
position  and  time.  Specifically,  the  stick  distance  (Figure  39)  is 
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Thus,  we  believe  the  most  practical  approach  will  be  a fixed  stick 
spacing  together  with  minor  (time-varying)  point-^head/behiud.  We  have 
examined  the  residual  isoplanatic  degradation  resulting  from  choosing 
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Xstick  at  comPromise  9 angles  of  22.5°  and  30°,  at  A = 3 pm,  and  found 
it  to  be  small.  It  will  be  more  serious  at  visible  wavelengths,  but  we 
have  yet  to  quantify  it. 

In  any  event,  given  a fixed  stick  length,  one  is  faced  with  the 
combined  problem  of  generating  angular  and  frequency  offsets  in  a system 
that  generally  does  not  accommodate  them  in  the  e’ementary  forms  illus- 
trated in  Figure  37.  However,  we  have  founu  that  it  is  possible  to 
generate  these  offsets  by  several  approaches  (as  discussed  below).  We 
call  these  subsystems  "override"  systems  since  they  override  the  natural 
tendency  of  the  conjugator. 

In  spite  of  the  Doppler  and  point-ahead  problems,  which  are  incurred 
as  a consequence  of  the  orbital  motions,  there  are  offsetting  benefits. 
Since  the  reference  oscillator  or  corner  is  well  ahead  of  the  relay 
mirror  by  the  time  the  high-power  phase-conjugated  beam  arrives  at  the 
relay  mirror  (Figure  39(d)),  the  feedback  coupling  into  the  oscillator 
(system  A)  or  the  reference  corner  (system  B)  is  minimized.  If  this 
angular  spatial  separation  is  not  adequate  to  suppress  the  main  beam 
feedback  to  the  desired  levels,  the  polarization  separation  techniques 
previously  discussed  in  the  OSC/MOPA  compensation  sections  can  be  used. 
Note,  however,  that  this  decoupling  presents  a problem  to  the  cw  iterated 
pulse  system  (system  C of  Figure  37).  In  particular,  we  desire  to  form 
a well-collimated,  mildly  converging  beam  centered  on  the  relay  mirror, 
yet  we  must  simultaneously  illuminate  the  reference  with  sufficient 
energy  to  return  a pulse  of  adequate  strength  to  the  ground  system  to 
establish  the  next  iterated  pulse.  A rough  calculation  suggests  that 
there  is  not  adequate  energy  in  the  sidelobe  structure  of  the  main  beam 
to  give  the  feedback  desired  for  the  system  (typically  we  would  design 
for  2%  total  feedback  with  10%  of  the  main  beam  being  extracted  for  Lhe 
reference  retro).  We  have  identified  two  possible  solutions  to  this 
problem.  One  approach  is  to  extract  a portion  of  the  beam  from  the 
vicinity  of  the  relay  mirror  and  to  transfer  it  to  the  reference  site 
by  a mini  relay  system.  The  other  approach  is  to  split  the  override 
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Figure  39.  Time  progression  of  a pulsed  wavefront  in  a carrot-on-a- 

stick  retro  system.  In  the  ideal  system  the  initial  pulse 
is  launched  to  point  ahead  of  the  corner  reference  (a), 
intercepts  it  and  is  retro  reflected.  (b),  is  reflected 
and  conjugated  at  the  ground  station  (c),  and  intercepts 
the  relay  mirror  on  axis  (d).  Since  the  conjugated  wave- 
front  retraces  exactly  the  same  path  the  isoplanatic  prob- 
lem is  minimized. 


system  in  such  a way  that  it  generates  two  beams,  one  (90%  of  the 
energy)  centering  on  the  relay  and  the  other  (10%  of  the  energy)  center- 
ing and  focusing  on  the  retro.  Either  approach  extracts  the  reference 
with  high  efficiency.  With  the  simpler  split  beam  approach,  the  beam 
centered  on  the  retro  will  be  pointed  ahead  and  thus  will  traverse  a 
somewhat  different  atmospheric  path  (i.e.,  it  will  be  subject  to 
isoplanatic  effects).  This  does  not  matter  for  small  retros. 

Systems  B and  C can  obtain  improved  oscillator  injection  efficiency 
by  means  of  a polarization-sensitive  beam  splitter  and  a polarization 
rotator  on  the  satellite.  For  this  purpose,  a quarter-wave  plate  backed 
by  a cats-eye  retro  system  could  replace  the  corner  reflector.  In  this 
application,  the  phase  conjugator  should  be  designed  to  reflect  with  the 
inci.  ' polarization  state. 

C.  FOCUS,  TILT,  AND  DOPPLER  OVERRIDE  SYSTEMS 

As  discussed  in  the  previous  section,  the  natural  tendency  of  a 
basic  nonlinear  phase  cpnjugator  is  to  return  the  energy  at  the  effective 
frequency  of  the  reference  to  an  aimpoint  centered  on  the  reference,  with 
a focus  corresponding  to  the  divergence  exhibited  in  the  radiated  wave- 
front  of  the  reference.  For  those  cases  where  the  reference  is  moving 
and  is  not  located  on  the  target,  the  basic  system  may  miss  the  target 
as  well  as  one  or  more  additional  key  system  objectives.  More 
specif ically : 

• The  motion  of  the  target  or  of  a relay  satellite  in  a relay 
system  during  the  reference  pulse  transit  time  may  provide 
a complete  miss. 

• For  systems  like  the  relay  systems  where  the  reference  is  not 
on  the  target,  the  natural  focus  provided  by  the  conjugator 
is  generally  incorrect  for  focusing  on  the  target.  Even  with 
target  referencing,  the  target  motion  during  the  reference 
propagation  time  may  produce  some  focusing  error. 

• The  Doppler  shift  of  the  reference  or  the  target  may  cause 
the  returning  wavefront  and  the  conjugated  wavefront  to  be 
sufficiently  shifted  in  frequency  to  fall  outside  the  gain 
profile  of  the  lasing  medium. 
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Since  these  problems  would  significantly  degrade  the  performance  of 
our  projected  system,  we  believe  that  some  type  of  "override"  on  the 
basic  conjugator  operation  will  be  required  to  compensate  for  them.  We 
have  identified  two  basic  approaches  to  such  override  systems.  First, 
the  outgoing  and  return  paths  can  be  effectively  separated  such  that 
tilt,  focus,  and  frequency  changes  can  be  introduced  on  only  the  outgoing 
path.*  As  discussed  below,  the  path  separation  may  be  either  temporal 
(limited  to  short  pulse  systems)  or  spatial  (via  polarization  separation). 
Second,  an  appropriate  perturbation  can  be  introduced  in  one  or  more 
pump  angles  and  frequencies  such  that  the  retro  wavefront  retains  its 
conjugated  properties  but  has  a shifted  frequency,  tilt,  or  focus  as  a 
consequence  of  the  pump  perturbations. 

1 . Override  by  Path  Separation 

This  subsection  discusses  temporal  and  spatial  techniques  for  focus 
and  tilt  override.  Frequency  override  can  in  principle  also  be 
incorporated  by  including  a frequency-shifting  unit  (such  as  an  accounts 
optic  cell),  but  this  will  not  be  discussed  in  detail. 

Temporal  techniques  require  short  pulse  operation,  and  we  will  use 
the  three-pulse  satellite  system  discussed  above  as  an  example.  We  have 
added  a gated  override  box  to  the  system  in  Figure  40(a),  and  in  Fig- 
ure 40(b)  we  have  illustrated  one  approach  to  constructing  the  gated 
override.  Basically,  just  following  the  passage  of  the  tail  end  of  the 
return  probe  pulse  through  the  electrooptical  element,  a high-speed 
switch  (such  as  a thyratron)  transfers  a fixed  voltage  (and  charge)  onto 
the  electrooptical  electrodes  via  a capacitor  precharged  to  a voltage 
appropriate  to  a prescribed  tilt  and  focus  override.  The  tilt  system 
consists  of  cascades  of  prisms  with  alternating  directions  of  the  trans- 
verse electric  field  on  successive  prisms.  A one-dimensional  focusing 
exement  that  uses  the  alternation  principle  is  also  illustrated.  The 


* 

The  outgoing  path  is  conceptually  the  more  natural  location  for  the 
override  system,  but  the  incoming  path  will  serve  equally  well  with 
inverted  (reversed)  override  perturbations  to  those  actually  required. 
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extension  of  the  focus  system  to  two  dimensions  involves  some  complexity, 
which  takes  us  beyond  our  present  scope,  but  solutions  have  been  worked 
out.*  The  electrodes  of  the  system  are  discharged  following  the  transit 
of  the  third,  or  high-energy,  excitation  pulse.  The  discharged  process 
is  allowed  to  be  more  leisurely  since  it  need  only  be  completed  before 
the  onset  of  the  next  oscillator  initiation  pulse  of  the  next  three- 
pulse  sequence.  If  the  corner  movement  does  not  quench  the  potential 
fourth  pulse,  the  deflection  potential  of  the  electrooptical  element, 
possibly  with  an  associated  pinhole,  can  be  used  to  quench  the  fourth 
pulse  return  to  the  conjugator.  For  example,  this  could  drive  it  below 
the  threshold  of  an  SBS  conjugator. 

In  the  second  class  of  path-separation  systems,  illustrated  in 
Figure  41,  a four-wave-mixer  conjugator  is  used.  It  switchees  the 
(linear)  polarization  state  of  the  incident  signal  wavefront.  This 
permits  the  outgoing  and  returning  wavefronts  to  be  effectively  separated 
with  high  efficiency.  The  polarization  state  reflected  off  of  the 
reference  retro  system  is  similarly  switched  in  this  instance  by  a 
double-pass  through  a quarter-wave  plate.  Thus,  the  outgoing  wavefront 
is  everywhere  horizontally  polarized,  and  the  return  wavefront  is  every- 
where vertically  polarized.  The  brewster-angle  multiplayer  filter  is 
only  one  of  several  possible  techniques  for  achieving  polarization 
separation.  At  3 pm  or  longer  wavelengths,  diffraction  grating  systems 
will  undoubtedly  ^e  the  preferred  approach. 

2 . Pump  Perturbation  Override  Systems 

We  have  found  that  it  is  also  possible  to  perturb  the  wavelengths 
And  angles  in  the  pump  source  of  a four-wave  mixing  conjugator  in  such 
a way  that:  (1)  the  Doppler  shifts  are  cancelled  in  the  phase-conjugated 
wavefront  by  the  pump  frequency  perturbations,  (2)  the  conjugated  wave- 
front  is  returned  with  the  desired  small  angular  offset  to  yield  the 


*Hc 

For  example,  see  J.F.  Lotspeich,  U..  Patent  No.  3,892,469,  "Electro- 
optical  Variable  Focal  Length  Lens  Using  Optical  Ring  Polarizers," 
July  1,  1975. 
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Figure  41.  Focus  override  by  path  splitting  by  polarization 
rotation  via  pump  manipulation. 
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point-ahead  adjustment,  (3)  exact  k-vcctor  matching  (for  efficient 
mixing)  :s  retained,  and  (4)  at  least  a partial  focus  override  correction 
can  be  incorporated.  The  system  configurations  that  permit  these  over- 
rides of  the  basic  conjugator  operation  are  discussed  in  Appendix  D. 

In  addition,  the  sizes  of  the  required  pump  perturbations  are  quantified 
in  the  appendix.  The  approach  will  accommodate  either  a full  point- 
ahead  correction  or  a part'al  correction  with  the  carrot-on-the-st ick 
approach . 

D.  I SOPl.ANATIC  EFFECTS  ASSOCIATED  WITH  THE  EXTENT  OF  THE  REFERENCE 

As  mentioned  in  the  introduction,  even  when  the  reference  is  moved 
ahead  of  the  relay  by  exactly  the  point-ahead  angle,  there  are  additional 
isoplanatic  problems  that  occur  with  a relay  system  operating  in  law 
orbits.  These  problems  are  associated  with  the  fact  that  it  is  generally 
impractical  to  make  the  reference  with  an  extent  comparable  to  that  of 
the  relay  mirror.  In  particular,  as  illustrated  in  Figure  42,  it  is 
advantageous  to  use  a reference  that  has  a small  aperture  to  hold  down 
system  cost  and  to  minimize  OPD  errors  associated  with  the  reference 
optics  (which  would  be  impressed  on  the  main  laser  output  with  any 
type  of  conjugation  system).  As  illustrated,  from  a ray  optics  viewpoint, 
there  is  an  angular  shift  Op  between  the  reference  rays  and  the  beam 
rays  that  is  typically  10  to  20  prad  at  the  outer  rays  for  low-altitude 
orbits.  A broad  comparison  of  these  numbers  to  the  results  given  in 
Figure  D-4  of  Appendix  1)  strongly  suggests  that  this  class  of  isoplanatic 
effect,  which  we  call  "focus  mismatch,"  will  be  troublesome  at  visible 
wavelengths.  There  does  not  appear  to  be  a significant  focus  mismatch 
problem,  even  for  point  source  references,  for  the  synchronous  orbit 
problem  since  the  offset  angles  are  too  small. 

Although  we  believe  this  is  an  important  problem  deserving  of  further 
stuay,  it  appears  to  be  a basic  problem  associated  with  referencing 
from  relay  systems  at  visible  wavelengths  and  one  that  is  independent 
of  the  adaptive  optics  system  employed,  be  it  an  NPC  or  a conventional 
svstem. 
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We  have  obtained  quantitative  measurements  on  the  correction  of  severely  aberrated  laser  1 jams  using  stimulated 
Brillouin  scattering  (SBS)  ut  0.69  «im.  We  have  shown  that  under  certain  conditions  SBS  can  be  used  to  restore 
an  aberrated  optical  beam  to  ita  original  unaberrated  condition.  When  an  optical  beam  double  passes  an  aberrat- 
ing region  after  reflecting  from  an  “ordinary”  mirror  (i.e„a  plane  mirror)  the  aberration  is  twice  that  obtained  from 
a single  pass.  However,  when  the  aberrated  beam  enters  a medium  that  allow  SBS  to  occur,  it  emerges  from  ita 
second  pais  through  the  aberrating  medium  In  the  same  condition  as  that  in  which  it  originally  entered.  Quantita- 
tive experiments  are  described  in  which  a single-mode  ruby  laser  beam  is  Intentionally  aberrated  by  passing  it 
through  an  etched  plate.  When  the  beam  la  allowed  to  double- pees  the  plate  using  an  ordinary  reflector  (i.e.,  plane 
mirror),  the  beam  divergence  is  more  than  10  tlmaa  the  diffraction-limited  divergence.  However,  when  we  replace 
the  ordinary  reflector  with  a call  in  which  SBS  can  take  place,  the  SBS  reflected  beam  is  restored  to  diffraction-lim- 
ited diver^i.ice  when  it  la  allowed  to  pass  back  through  the  aberrating  medium.  Application*  of  this  time-reversal 
or  phase-reversal  technique  for  correcting  aberrations  in  optical  trains  and  atmospheric  turbulence  are  discussed. 


The  correction  of  aberrations  introduced  by  the  atmo- 
sphere or  by  optical  components  in  the  propagation  of 
laser  beams  is  of  great  interest  and  current  activity.1*2 
Adaptive  optics  techniques  currently  under  develop- 
ment for  reducing  the  effects  of  atmospheric  turbulence 
and/or  optical  train  distortion!  on  laser  beams  [e.g., 
coherent  optical  adaptive  techniques  (COAT)]  involve 
systems  that  adjust  the  phasefront  of  the  transmitted 
beam  to  compensate  for  p.  ie  aberration  by  monitoring 
backscatter  from  the  target.  The  correction  is  accom- 
plished by  ui  ing  some  form  of  discrete,  multichannel, 
phasefront  corrector  driven  by  electronic  servos,  Buch 
as  a deformable  mirror. 

In  contrast,  certain  nonlinear  optical  interactions  can 
generate  the  spatial  phase  conjugate  of  a distorted 
wavefront,  which  then  can  be  transmitted  through  the 
original  distorting  optical  path  to  form  a corrected  beam 
(the  return  beam  from  a glint  upon  a target  contains  all 
the  necessary  information  to  correct  for  the  atmo- 
sphere’s distortion).  These  nonlinear  interactions 
automatically  perform  a phasefront  correction  that  is 
spatially  continuous  over  the  entire  cross  section  of  the 
beam  without  any  external  wavefront  sensing  or  elec- 
tronic controls  and  may  well  have  an  economy,  Bpeed, 
and  simplicity  beyond  conventional  adaptive  optical 
schemes. 

The  work  described  here  represents  the  first  report 
of  a quantitative  experimental  evaluation  of  phase 
conjugation  derived  through  a nonlinear  optical  effect 
We  also  describe  a concept  where  this  effect  can  be  used 
to  correct  the  aberrations  of  a laser  transmitter  props  • 
gated  through  a turbulent  atmosphere. 

The  basic  concept  of  nonlinear  phase  conjugation  for 
laser  beam  correction  through  the  atmosphere  is  shown 
schematically  in  Fig.  1.  The  process  may  be  described 
as  occurring  in  the  following  steps: 


Fig.  1.  A conceptual  pulsed  laser  system  using  SBS  for  cor- 
rection of  atmospheric  and  laser-induced  distortions. 


1.  The  first  transmitted  pulse  propagates  to  the  tar- 
get, assumed  for  the  purpose  of  this  example  to  be  a 
Bingle  unresolved  glint  The  purpose  of  this  first  pulse 
is  to  illuminate  the  glint  target,  which  then  serves  as  a 
test  source  at  the  wavelength  of  interest. 

2.  The  light  reflected  from  the  glint  propagates 
through  the  distorting  medium  toward  the  transceiver, 
arriving  as  an  aberrated  wavefront. 

3.  The  phase  conjugate  of  this  distorted  return  wave 
is  generated  by  a nonlinear  optical  conjugator. 

4.  After  coherent  amplification  to  a desired  power 
level,  the  phase-corrected  pulse  retraverses  the  dis- 
torting atmosphere,  which  now  restores  its  phase  co- 
herence so  that  the  entire  beam  is  focused. 

The  weak  backscatter  to  be  expected  from  a distant 
target  would  suggest  that  a very-high-gain  amplifier 
would  be  necessary  to  overcome  the  losses.  In  the  last 
few  years  the  technology  for  this  high  gain  has  been 
developed  for  1.06-  and  10.6-^m  laser-fusion  systems. 
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A gain  of  nearly  60  dB  has  been  demonstrated,  and  this 
appears  consistent  with  operation  over  a range  of  several 
kilometers.3 

An  example  of  a nonlinear  process  that  possesses  the 
ability  to  create  directly  the  wave  whose  Held  is  the 
complex  conjugate  of  a signal  wave  is  stimulated  Bril- 
louin  scattering  (SBS).  It  was  first  noticed  in  19664  that 
the  backscatter  from  SBS  possessed  the  peculiar 
property  that,  if  the  incident  laser  beam  was  diverging, 
the  backscattered  light  formed  a converging  beam  of 
matching  cone  angle.  Zel'dovich  etal . 6 and  Nosach  et 
al.6  have  reported  the  ability  of  SBS  to  compensate  for 
aberrations  in  a ruby  laser  beam.  Zel’dovich6  and  re- 
cently several  other  Soviet  researchers7  have  given 
semiquantitative  arguments  that  the  backscattered 
light  of  SBS  can  everywhere  retrace  the  path  of  the  in- 
cident light  and  is  equivalent  to  the  spatial  phase  con- 
jugate wavefront  of  the  incident  light.  The  conjugate 
field  can  be  viewed  as  equivalent  to  the  incident  field 
traveling  backward  in  time.8,9  Nosach6  used  this  effect 
to  compensate  for  aberrations  in  a ruby  amplifier 
rod. 

In  view  of  the  importance  this  process  migT‘  have  to 
adaptive  optics  and  because  of  theoretical  uncertainties 
accompanying  our  present  understanding  of  this  pro- 
cess, we  have  prepared  the  experiment  shown  in  Fig.  2, 
following  closely  the  work  of  Nosach  et  a/.,6  especially 
in  the  design  of  the  highly  multimode  waveguide  CS2 
cell.  A TEMoo  single  longitudinal-mode  ruby  oscillator 
passively  Q switched  by  a dye  cell  provides  an  output 
pulse  of  17-nsec  duration.  This  is  amplified  and  then 
controlled  by  an  adjustable  attenuator  to  deliver  about 
25  mJ  to  the  SBS  cell. 

An  aberration  can  be  purposely  introduced  into  the 
beam  at  position  2 or  3 in  Fig.  2 to  test  the  corrective 
properties  of  SBS  backscatter.  The  aberrator  is  a mi- 
croscope slide  etched  in  HF.  The  original  beam  0,'iality 
is  monitored  by  placing  a partially  reflecting  mirror  in 
position  2 and  reflecting  a portion  of  the  light  back  into 
the  1-m  focal -length  camera  through  the  beam  splitter. 
The  degree  of  aberration  is  measured  by  placing  the 
aberration  in  front  of  the  reference  reflector  at  position 
2,  causing  the  light  to  double-pass  the  aberrator.  The 
SBS-corrected  process  is  accomplished  by  passing  the 
light  through  the  amplifier  and  attenuator,  through  the 
aberrator  now  at  position  3,  and  into  the  focusing  lens 
and  CS2  waveguide.  The  lens  is  used  to  focus  the  beam 
to  a waist  just  in  front  on  the  cell  entrance  window  giv- 
ing a divergent  entrance  beam.  Since  the  index  of  CS2 
is  higher  than  the  glass,  this  diverging  beam  is  com- 
pletely internally  reflected  in  the  2.5-mm-i.d.  1-m-long 
capillary.  The  waveguide  SBS  cell  serves  to  increase 
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Fig.  2.  Experiment  to  measure  the  beam  divergence  of  SBS 
backscatter  at  0.69  >im. 
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Fig.  3.  Measured  full-angle  beam  divergences  of  original, 
aberrated,  and  corrected  beams  and  their  appearance  in  the 
far  field. 


the  interaction  length  and  percentage  of  light  hack- 
scattered  without  increasing  power  density  to  the  point 
where  1 ir-induced  breakdown  of  CS2  can  occur  and 
without  increasing  the  total  power  input  to  the  point  of 
self-focusing.  The  slightly  redshifted  SBS  backscatter 
retraces  the  path  of  the  original  beam  through  the 
multiple  bounces  in  the  waveguide  and  through  the 
aberrator  and  the  amplifier.  A portion  is  sampled  by 
the  beam  splitter  and  the  1-m  camera. 

Divergence  of  these  beams  is  measured  by  a technique 
using  a series  of  exposures  at  differing  attenuations  at 
the  focal  point  of  a 1-m  lens.10  In  this  way  the  actual 
profile  of  the  original  aberrated  and  corrected  beams 
can  be  determined  independently  of  film  linearity.  In 
a plot  of  the  log  of  the  relative  intensity  versus  spot  di- 
ameter squared,  as  recorded  on  film  of  sharp  threshold, 
a Gaussian  beam  will  result  in  a straight  line  plot,  with 
slope  proportional  to  the  1/e  Gaussian -beam  diame- 
ter. 

Beam  profile  plots  were  taken  for  the  reference,  ab- 
errated, and  SBS-corrected  beams,  and  all  were  ap- 
proximately Gaussian.  The  full-angle  beam  diver- 
gences obtained  from  the  slope  of  these  plots  are  0.57 
± 0.06, 6.6  ± 0.9,  and  0.44  ± 0.06  mrad.  The  SBS  pro- 
cess provides  a corrected  beam  that  has  a divergence 
almost  indistinguishable  from  the  original  beam  and 
substantiates  earlier  results.6,6  In  fact,  the  SBS  beam 
divergence,  0.44  mrad,  is  the  calculated  diffraction- 
limited  divergence  of  this  laser  beam.  (The  slight  im- 
provement over  the  reference  beam  is  unexplained  but 
could  result  from  the  reference  retroreflector's  being  out 
of  flat  by  X/6  across  the  beam.)  In  Fig.  3 t he  normalized 
beam  intensity  is  plotted  versus  divergence  or  radius  at 
thd  focal  plane  of  the  camera.  Also  shown  are  photo- 
graphs of  the  reference,  aberrated,  and  corrected  beams. 
The  corrected  beam  has  horizontal  sidelobes  not  present 
in  the  reference.  The  presence  of  these  sidelobes  makes 
a negligible  contribution  to  the  overall  beam  divergence 
since  their  intensity  is  very  low  (the  outermost  two  data 
points  in  the  plot  of  the  corrected  beam  divergence  are 
the  result  of  these  sidelobes). 

There  is  some  theoretical  justification  for  concern 
that  other  wavefronts  in  addition  to  the  desired  phase 
conjugate  wavefront  may  exist  in  the  backscattered 
beam.  For  example,  a plane  wavefront  backscattered 
by  the  SBS  might  go  undetected  after  passing  through 
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the  aberrator  if  the  resulta  it  diffused  beam  gave  a 
maximum  intensity  below  the  film’s  threshold  of  sen- 
sitivity. Locating  the  source  of  aberration  ahead  of  the 
beam  splitter  and  adjacent  to  the  oscillator  did  not  re- 
sult in  a bright  focused  spot  at  the  camera  focal  plane, 
which  would  be  the  result  of  a collimated  return  wave. 
We  are  now  in  the  process  of  exploring  the  effect  of 
power  and  interaction  length  upon  the  efficiency  of  the 
correction. 

We  have  also  found  that  an  optical  waveguide  appears 
necessary  to  the  phase  conjugation  process.  When  the 
beam  was  focused  within  a conventional  cell  with 
comparable  total  input  power  and  comparable  per- 
centage backsi  atter,  the  backscattered  light  was  not 
capable  of  correcting  the  divergence  introduced  by  our 
aberrator  but  ret  trned  a reflection  with  divergence 
similar  to  an  ordinary  reflection.  Thus  the  selection 
process  that  allows  ascendency  of  the  phase  conjugate 
phasefront  over  all  other  possible  phasefronts  seems  to 
be  dependent  on  propagation  of  confined  modes  within 
a waveguide. 

An  oversimplified  but  intuitively  satisfying  physical 
picture  of  this  process  is  illustrated  in  Fig.  4.  The  phase 
aberrated  beam  is  represented  by  a wavefront  with  a 
simple  step.  An  ordinary  reflection  results  in  the 
doubly  aberrated  wavefront  shown  as  Fig.  4.  As  a result 
of  the  SBS  interaction,  the  aberrated  pump  wave 
creates  aberrated  hypersonic  waves  with  identical 
phasef  ronts  (as  in  Fig.  4,  bottom)  which  act  as  a moving 
"deformable”  dielectric  reflector  that  yields  the  con- 
jugate s -attered  wave. 

Other  phase  conjugating  processes  that  have  been 
proposed  include  parametric  down-conversion11  and 
four-wave  mixing.912  At  this  time  SBS  appears  to  be 
quite  efficient,  yielding  from  10  to  70%  backscatter, 
depending  on  the  degree  of  aberration  introduced. 

Other  promising  applications  of  this  class  of  correc- 
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Fig.  4.  A simplified  view  of  the  SBS  phase  conjugating  pro- 
cess as  a moving  multilayer  reflector. 


tive  techniques  include  correction  of  large-diaineter 
laser  resonator  cavities  as  well  as  large  optical  trains, 
such  as  those  proposed  for  laser  fusion  or  space  op- 
tics. 

In  summary,  we  have  presented  the  first  report  of  a 
quantitative  evaluation  of  SBS  phase  conjugation  and 
suggested  a means  by  which  phase  conjugation  can  be 
used  for  adaptive  optics. 

We  wish  to  acknowledge  fruitful  discussions  with  V. 
Evtuhov,  R.  W.  Hellwarth,  D.  M.  Pepper,  and  A.  Yariv 
and  thank  T.  Horne  for  his  technical  assistance. 
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We  consider  the  backward  optical  wav,-  stimulated  by  a multimode,  monochromatic,  incident 
optical  wave  in  a waveguide  filled  with  a transparent  nonlinear  medium,  when  the  incident  wave 
is  negligibly  perturbed  by  the  nonlinear  proceeaea.  We  derive  the  conditions  on  guide  length,  area, 
mode  number,  and  Stokes  shift  in  order  that  a given  high  percentage  of  the  power  in  the  backacat- 
tered  field  be  the  "phase  conjugate”  of  the  incident  field,  i.e.,  be  proportional  to  its  complex  conjugate 
in  the  entrance  plane  of  the  waveguide. 


I.  INTRODUCTION 

When  a strong  monochromatic  wave  (at  r)  is  incident  on 
a transparent  medium  it  causes  waves  at  lower  frequencies  u 
to  experience  exponential  gain  if  their  frequency  offset  cor- 
responds to  the  frequency  of  some  excitation  in  the  medium. 
If  this  excitation  is  an  acoustic  wave,  the  effect  is  called 
stimulated  Brillouin  scattering  (SBS),  otherwise  it  is  called 
stimulated  Raman  scattering  (SRS) . It  has  long  been  known 
that  laser  sources  may  produce  stimulated  gain  that  is  large 
enough  (~1012)  so  that  spontaneously  scattered  light  experi- 
ences sufficient  amplification,  in  a single  pass  out  of  the  in- 
cident-beam region,  toemerge  with  a power  approaching  that 
of  the  incident  beam.  More  recently,  it  has  been  observed 
that,  when  the  incident  beam  is  multimode  (i.e.,  has  a complex 
wave  front),  the  high  intensity  backscattered  SBS  or  SRS 
waves  generated  in  this  single- pass  process  may  be  nearly  the 
‘‘phase  conjugate”  to  the  incident  wave.1-6  That  is,  in  the 
entrance  plane,  the  complex  amplitude  E*  (r)  of  the  stimu- 
lated wave  is  nearly  equal  to  the  complex  conjugate  E)  (r)  of 
the  incident  wave  times  a constant.  When  the  incident  and 
scattered  frequencies,  v and  <u,  are  nearly  equal,  phase  con- 
jugation makes  the  scattered  wave  appear  to  be  nearly  a 
time-reversed  replica  of  the  incident  wave  in  a large  region  of 
space.  Zel'dovich  et  al. 1 suggested  that,  when  w ~ v,  this  ef- 
fect should  come  from  the  solution  of  the  nonlinear  Maxwell 
equations  after  they  were  linearized  with  respect  to  the  scat- 
tered fields.  They  pointed  out  that  the  resulting  complicated 
set  of  coupled  linear  equations  for  the  spatial  amplitudes 
would  have  various  solutions  which  exhibited  spatial  gain,  and 
that  one  solution  should:  (a)  have  F „ ~ E’  in  some  plane,  and 
(b)  have  significantly  (“2  to  3 timet”)  higher  gain  than  the 
others.  As  the  SBS  phase-conjugavion  experiments  had  been 


performed  in  an  optical  waveguide,  Sidorovich3  examined  the 
equations  appropriate  to  an  optical  waveguide,  and  agreed 
with  these  conclusions,  however,  surmising  some  “necessary” 
conditions  for  this  effect  for  which  we  will  show  counter  ex- 
amples here.  Zel'dovich  and  Shkunov  pointed  out  that,  when 
the  beams  interact  in  free  space  (rather  than  in  a waveguide), 
phase  conjugation  can  occur  even  when  the  frequency  sepa- 
ration between  the  beams  was  large,  as  in  SRS.4  However, 
the  conditions  they  find  differ  from  those  we  will  find  here. 
Subsequently,  Zel’dovich  et  al*  observed  SRS  from  the  656 
cm-1  Raman  line  of  CS2  to  be  largely  a phase  conjugate  to  an 
incident  beam  comprising  about  300  transverse  modes.  Wang 
and  Giuliano  verified  the  high  degree  of  phase  conjugation 
present  in  SBS.6 

In  this  paper  we  will  analyze  the  degree  of  phase  conjugation 
present  in  waves  which  are  stimulated  (from  noise)  in  a 
waveguide  by  a strong  multimode  beam,  whose  propagation 
is  assumed  to  be  negligibly  disturbed  by  the  stimulated  pro- 
cesses and  negligibly  attenuated  by  linear  losses.  Phase 
conjugation  by  stimulated  scattering  from  an  unguided  beam 
(which  requires  more  power  and  may  be  spoiled  by  self-fo- 
cusing, breakdown,  etc.)  may  be  approximated  by  considering 
a waveguide  whose  length  equals  the  diffraction  length  of  the 
incident  beam.  We  find  that,  even  when  a gap  exists  between 
the  gain  of  one  solution,  or  wave  pattern,  and  that  of  the 
others,  the  maximum-gain  solution  is  never  a perfect  phase 
.conjugate  to  the  pump  (input)  wave  pattern.  However,  a very 
large  fraction  (>90%)  of  this  maximum-gain  solution  is  phase 
conjugate  to  the  pump  under  conditions  which  are  much  less 
restrictive  than  those  suggested  by  previous  workers.3  4 First, 
the  waveguide  must  have  a finite  number  of  modes  that  can 
be  excited  and  the  transverse  field  patterns  of  the  pump  and 


1050  J.  Opt.  Soc.  Am.,  Vol.  68,  No.  8,  August  1978  0030-3941/78/6808-1050*00.50  © 1978  Optical  Society  of  America  1050 


(3) 


stimulated  modes  must  be  essentially  the  same.  This  con- 
stitutes very  little  restriction  in  practice.  More  importantly, 
the  product  of  the  waveguide  length,  the  Stokes  wavelength 
shift,  and  the  number  of  excited  waveguide  modes  divided  by 
their  cross-sectional  area  must  be  less  than  a number  that  is 
of  the  order  of  the  square  root  of  the  maximum  acceptable 
fraction  that  is  nut  phase  conjugate. 

Unlike  Ref.  3,  we  find  (a)  that  large  differences  in  the  am- 
plitudes of  the  excited  incident  modes  do  not  spoil  phase 
conjugation;  lb)  the  waveguide  does  nut  have  to  be  much 
longer  than  a diffraction  length;  (c)  the  incident  beam  solid 
angle  dees  not  have  to  exceed  the  nonlinear  index  change. 
The  free-space  (no  guide)  theory  of  Ref.  4 also  produces 
somewhat  different  conditions  that  depend  on  the  value  of  the 
stimulated  gain. 

Unlike  previous  analyses,  our  conclusions  are  based  on  a 
general  perturbation  theory  for  the  non-phase-conjugate 
fraction  in  the  wave  of  maximum  gain,  valid  for  any  distri- 
bution of  incident-wave  amplitudes  and  phases.  We  apply 
this  theory  to  several  models  in  a rectangular  waveguide  which 
have  arbitrary  distribution  of  phase  and  arbitrary  distribution 
among  four  possible  arbitrary  amplitudes.  What  seems  truly 
remarkable  is  that  the  many  extraneous  terms  in  the  original 
coupled-mode  equations  fail  to  spoil  the  phase-conjugate 
nature  of  the  guided  wave  with  maximum  gain,  at  least  in  the 
quite  representative  classes  of  input  waves  in  a rectangular 
waveguide  which  we  considered.  A general  proof  of  the  con- 
ditions for  phase  conjugation  of  all  possible  inputs  into  all 
possible  guide  shapes  is  still  lacking.  However,  we  feel  that 
our  conditions  are  accurate  for  cases  of  interest  in  the  labo- 
' ratory. 

N.  FORM  OF  THE  NONLINEAR  POLARIZATION 

We  will  consider  the  interaction  in  a straight  waveguide  of 
a multimode  “pump”  wave  or  input  wave  with  stimulated 
baclucattered  “Brillouin"  waves,  i.e.,  waves  generated  by  in- 
teractions with  backward  acoustic  waves.  The  pump  waves 
will  be  assumed  to  be  monochromatic  at  angular  frequency 
» and  to  be  Fourier  analyzable  inside  the  guide  in  terms  of 
plane  waves  whose  wave  vectors  lie  mainly  in  a small  cone  of 
half-angle  0 about  the  waveguide  axis  z.  We  will  look  for  so- 
lutions of  the  nonlinear  Maxweli  equations  which  are  back- 
ward-scattered waves  of  a single  (Brillouin-shifted)  frequency 
<■>;  solutions  at  different  frequencies  superpose.  However, 
stimulated  gain  will  make  those  solutions  predominate  which 
lie  in  a narrow  range  about  frequency 

u>  ■ r-Wfl,  (1) 

where 

wfi  a 2nv,v!c  (2) 

is  the  acoustic  frequency.  Here  v,  is  the  sound  velocity,  n is 
the  (linear)  refractive  index,  and  c is  the  velocity  of  light  in 
vacuum. 

We  will  find  that  the  backscattered  waves  will  have  wave 
vectors  lying  mainly  inside  a cone  whose  half-angle  is  also  0, 
the  angle  containing  the  pump  waves.  (In  fact,  the  back- 
scattered  waves  will  often  be  nearly  the  time-reversed  replica, 
or  “phase  conjugate,”  of  the  pump  wave.)  Here  we  will  con- 
sider only  the  case 
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where  Q - wh/^h  is  the  “Q  value"  of  the  backward  acoustic 
waves  and  is  the  linewidth  of  the  resonant  acoustic  re- 
sponse (and  of  the  spontaneous  backward  Brillouin-scattered 
light).  When  Eq.  (3)  holds,  it  is  readily  verified  that  the 
magnitude?  of  the  sound  wave  vectors,  excited  hy  pump  and 
scattered  light  waves,  lie  within  the  bandwidth  Awft/u,  in 
which  the  acoustic  response  is  nearly  constant— the  same  as 
at  exact  resonance.  The  independence  of  Brillouin  response 
to  wave  vector  in  this  case  implies  that  the  amplitude  at  w of 
the  nonlinear  optical  polarization  density  P"1  (r)  depends  only 
on  the  optical  field  at  the  same  position  r.  That  is,  we  can 
neglect  spatial  dispersion  when  Eq.  (3)  holds.  It  is  well  known 
that  on  resonance,  the  nonlinear  Brillouin  polarization  at  r 
has  the  form 

P£*  ■ -iG  E,  E|  • Ew,  (4) 

where  E,  and  Ew  are  the  (complex)  amplitudes  of  the  pump 
and  backscattered  waves.  On  resonance,  G is  a real  constant 
which  we  will  relate  later  to  the  usual  plane-wave  Brillouin 
gain  coefficient.  [When  w deviates  from  the  resonance  con- 
dition (1),  G becomes  complex  and  of  smaller  magnitude.) 
Other  nonlinear  effects,  such  as  SRS,  give  a third-order  po- 
larization density  that  is  largely  imaginary  and  spatially  local 
as  in  Eq.  (4),  and  our  theory  here  will  apply  to  these  also,  after 
possible  minor  modification  of  the  tensorial  character  of  Eq. 
(4). 

Although  Eq.  (3)  has  been  well-satisfied  in  all  experiments 
to  date,  we  note  that  it  need  not  be  if  the  waveguide-accep- 
tance angle  is  larger  than  Q~  m.  In  this  case  phase-conjuga- 
tion properties  can  be  predicted  by  considering  subsets  of 
excited  incident  modes,  each  of  which  interacts  with  a dif- 
ferent subset  of  phonon  waves,  the  results  of  each  subset  of 
interactions  being  treated  essentially  as  we  will  treat  the  case 
obeying  Eq.  (3). 

III.  FORMULATION  IN  A CYLINDRICAL 
WAVEGUIDE 

We  will  assume  that  the  nonlinear  polarization  density  does 
not  disturb  the  incident  pump  waves  so  that  their  complex 
amplitude  E„  (r)  can  be  written  inside  the  waveguide  as 

Ame„(xj)eik”"*.  (6) 

m 

Here  the  Am  are  complex  amplitudes  for  the  pump  wave  to 
be  in  various  normal  modes  of  the  guide,  whose  transverse 
mode  patterns  im(x,y)  are  normalized  so  that 

J*  dxdyi'm-in  * 6m„.  (6) 

We  will  assume  that  the  transverse  refractive  index  variations 
or  metallic  walls,  etc.,  forming  the  guide  are  such  that  the 
transverse  functions  im  do  not  differ  significantly  for  the 
various  temporal  frequencies  (pump  and  Stokes-shifted)  that 
are  propagating  in  the  guide.  However,  the  propagation 
constants  km,  are  frequency-dependent  and  obey 

ki,+ui-ny/c\  (7) 

where  u£  is  the  (frequency-independent)  eigenvalue  associ- 
ated with  the  function  im.  Here,  n„  the  refractive  index  at 
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the  guide  axis,  ia  auumed  to  be  teal,  i.e.,  linear  attenuation 
m the  interaction  regions  is  assumed  to  be  negligible. 

We  will  not  need  to  digress  here  on  the  characteristics  of 
waveguide  modes;  our  results  will  be  largely  independent  of 
the  detailed  nature  of  the  waveguide,  provided  that  it  prop- 
agates a number  of  modes  at » and  w whose  transverse  mode 
patterns  J are  congruent. 

IV.  NONLINEAR  MAXWELL  EQUATIONS  FOR 
MODE  AMPLITUDES 

We  now  consider  Maxwell's  equations  for  the  complex 
amplitude  E„.  of  the  waves  generated  in  the  guide,  at  frequency 
u.-.  under  the  influence  of  the  nonlinear  polarization  density 
(4).  We  will  try  a solution  of  the  nonlinear  Maxwell  equations 
of  the  form 

- Z Bn  (8) 

n 

which  is  similai  to  Eq.  (5)  except  that  the  waves  are  both 
traveling  and  growing  in  the  backwards  direction  (-z).  When 
there  is  no  nonlinear  term  (4)  and  y ■ 0,  Eq.  (8)  satisfies 
Maxwell's  equations  identically.  When  the  nonlinear  term 
(4)  is  added  as  a small  perturbation.  Maxwell's  equations  are 
satisfied  in  a length  L of  guide  provided  that 

yk  Bn  - 4*u >*c~2G  £ f f * dy  t ‘m . 4t  t)  • tn 

ijnJ  J 

X Kmijn  AjAj  Bn,  (9) 

where 

Kmijn  " *iA*f  dt/L  (10a) 

and  A*  * + k„  - kj,  - k„u.  It  is  useful  to  note  that 

|Kmi,ii|  "X_,siiu,  (10b) 

where  z s AkL/2  equals  zero  when  i ■ j and  m ■ n. 

in  Eq.  (9),  as  in  the  following  formulas,  repeated  space  in- 
dices are  assumed  to  be  summed  unless  otherwise  stated.  We 
have  approximated  a factor  hm„  on  the  left-hand  side  of  (9) 
by  its  average  k;  the  small  difference  is  negligible,  as  will  be 
seen  by  an  obvious  extension  of  the  perturbation  theory  which 
we  will  apply  to  more  important  terms  below. 

From  Eq.  (9)  it  is  seen  that  the  nonlinear  Maxwell  equations 
lead  to  a coupled  set  of  linear  equations  for  the  backward- 
scattered-wave  mode  amplitudes  B„,  when  these  are  too  small 
to  affect  the  incident-wave  amplitudes  Am.  We  have  not 
established  t.ie  properties  of  B„  for  all  possible  seta  of  Am,  but 
we  hav»  nolvea  for  B„  for  some  nontrivial  Am.  We  have  also 
derived  a perturbation  expression  for  the  non-phase-conju- 
gated  fraction  / of  that  solution  |B„|  of  Eq.  (9)  having  the 
largest  vulue  >o  of  gain  y.  We  also  calculate  ><>,  showing  that, 
when  / « 1,  >o  larger  enough  than  the  y for  any  other  solu- 
tion so  that  the  near-conjugate  solution  can  dominate  in 
practice.  The  perturbation  theory  needed  to  accomplish 
these  tasks  takes  a form  that  is  familiar  to  quantum  mechanics 
if  we  rewrite  Eq.  (9)  in  the  form  of  an  eigenvalue  equation, 
with  eigenvalue  Y,  in  a Hilbert  space  of  dimension  N equal 
to  the  number  of  transverse  mode  patterns  that  are  excitable 
in  the  guide.  This  reexpression  of  Eq.  (9)  is 


UHmn  + Vmn  + Umn)BnmYBm, 

ft 

where  the  “effective  unperturbed  Hamiltonian"  is 

HmnEa6mn  + fia"m  a„  (12) 

and  the  first  perturbation  Hamiltonian  is  (no  implied  sum- 
mations) 

Vmn  ■ Z (anl  “ <r)  0/  0/  &nm  + (0mn  “d)®*  arr  (13) 

I 

Here  the  a„  are  normalized  dimensionless  pump-beam  am- 
plitudes 

•*-4.  (DC  4.)  . (14) 

The  am„  and  0mn  matrices  are  derived  from  mode-overlap 
integrals  according  to 

amn  s SS  dxdy\e‘H‘*n\2  (15) 

and  (no  implied  summation) 

0mn  E SS  dxdy(6'm-in)*Km*mn.  (16) 

Their  averages  a and  0 are  defined  by 

® E L‘»mnl®m|2|®n|2  (17) 

mn 

and  similarly  for  0.  One  sees  that  the  a terms  in  Eqs.  (12)  and 
(13)  are  from  the  terms  in  (9)  for  which  i ■ j and  m ■ n;  the 
0 terms  are  from  the  terms  in  (9)  for  which  i ■ n and  j » m. 
Note  that  Hm„  and  Vmn  are  Hermitian. 

The  second  perturbation  Hamiltonian  is  (no  implied 
summation) 

Umn  E Z 0,a‘  J J dx  dy  - f,  f*  - 4nKnijn,.  (18) 

where  the  prime  on  the  summation  signifies  the  inclusion  of 
only  those  terms  not  included  in  Eqs.  (12)  and  (13).  The 
K„ujn  in  Eq.  ( 18)  have  Ak  0 except  when  degeneracies  occur 
which  make  either  At,-,  * kjr  and  k„v  * kmv  simultaneously, 
or  ki,  “ kn»  and  kJr  ■ kmu  simultaneously.7  In  practice,  the 
Kmij„  in  Eq.  (18)  for  which  Ak  a 0 can  be  omitted  from  con- 
sideration for  the  following  reasons.  First,  the  component  of 
E,  in  any  degenerate  manifold  of  guide  modes  can  be  taken 
as  one  of  the  basis  vectors  in  this  manifold.  Then,  no  terms 
will  exist  in  (18)  with  kif  * kjr.  The  second  case  is  possible 
when  wot  v,  uf  m u;,  and  uj  *u  jjj.  However,  then  it  is  usual 
that  e|i  • i j ■ 4)  • em  = 0.  In  a guide  with  high  symmetry,  the 
latter  may  not  bt!  the  case,  but  then  the  x-y  overiap  integral 
in  (18)  seems  to  be  greatly  reduced  over  that  in  (16)  arising 
from  i “ n and  jmm\ in  fact  it  reduces  to  zero  for  a rectangular 
waveguide.  Thus  we  will  neglect  the  terms  in  (18)  for  which 
Ak  is  accidently  near  zero,  considering  next  the  effects  of  the 
remaining  terms. 

If  the  difference  between  AkL  values  for  various  mode  in- 
dices is  larger  than  ir,  expressions  (10)  show  K„nj„  to  oscillate 
randomly  and  cause  a great  deal  of  cancellation  among  terms. 
If  the  area  of  the  guide  is  S,  and  there  are  N excitable  modes, 
then  u\  ~ N2/S,  « 1,  and  AkL  changes  by  at  least  t 
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for  most  changes  in  (f,  j,  n,m)\  hen 

NL  > Sk.  (19) 

We  expect  therefore  that  ( 'mn  can  be  neglected  compared  to 
Hmn  when  L > Ld/Nln,  where  i»  a diffraction  length 
{~-Sk/N  ,/2).  This  condition  is  met  by  free-spare  interactions 
and  all  practical  guides.  Therefore  we  proceed  to  solve  Eq. 
(ll)withUmn  ■!). 

Before  proceeding  it  is  useful  to  relate  the  eigenvalue  Y for 
any  N -dimensional  vector  solution  |W„|  of  Eq.  (11)  to  its  gain 
y and  the  commonly  defined  stimulated  gain  coefficient  n ( Np 
cm/MW)  for  plane  pump  and  scattered  waves.  The  “free- 
space"  coefficient  g is  the  same  as  that  lyS  + total  pump 
power  P)  in  a guide  that  is  excited  in  a single  mode  i?„  having 
a uniform  intensity  profile  («>,,  ■ const),  and  lor  which  Eqs. 
(6),  (15),  and  (16)  give  o = d * US.  In  this  case  V,,,,,  = Unm 
* 0 and  Bm  = a’n  is  the  solution  of  ( 1 1 ) with  eigenvalue 

Ypu.  - 2 i/S.  (20) 

Since  | A„  | 2 is  proportional  to  the  power  P in  the  guide,  a 
comnarison  of  Eqs.  (9)  and  ( 1 1)  shows  V <*  y/P  in  all  cases. 
Therefore  the  power  gain  y for  any  solution  of  (11)  with  ei- 
genvalue Y is  related  to  g b\ 

T(cm-1)  * Hi  Y(cm'2)g(cm/MW)  P(MW),  (21) 

where  the  most  common  units  used  are  indicated.  We  pro- 
ceed now  to  find  solutions  of  Eq.  (11)  and  their  gains  y for 
multimode  puinp  beams. 

V.  FIRST  APPROXIMATE  SOLUTION  FOR 
ARBITRARY  INPUT 

The  “unperturbed”  matrix  Hmn  in  (1 1)  is  of  the  particularly 
simple  form  of  a unit  matrix  plus  a projection  operator  Pa  onto 
the  conjugate  a’  of  that  vector  whose  components  a„  are  the 
complex  mode-amplitudes  of  the  incident  wave  in  the  guide. 
Therefore,  the  eigenstates  and  eigenvalues  G1'*  of  Hmn  are 
easily  seen  to  be  a‘,  with  eigenvalue 

G<n>  - a + d (22) 

and  any  set  of  vectors  b{n,]  orthogonal  toaj,  and  to  each  other, 
which  complete  the  space  of  guide  modes.  The  latter  have 
eigenvalues 

G<»>  - a.  (23) 

Note  that  when  all  guide  modes  have  the  same  linear  (or 
elliptical)  polarization,  and  when  u>  v,  so  that  by  Eq.  (10) 
Kmnmn  ~ * 1,  as  in  practical  SBS  experiments,  then 

5-?  (24) 

and 

G°  — 2G'.  (25) 

That  is,  the  conjugate  backscattered  wave,  whose  mode  am- 
plitudes are<C  would  experience  twice  the  gain  as  any  other 
mode,  provided  that  Vm„  can  be  neglected. 

We  will  find  that  whenever  the  scattered -wave’s  energy  is 
mostly  in  the  phase  conjugate  to  the  incident  wave,  its  gain 
is  nearly  twice  that  of  any  other  wave.  Furthermore  this 
condition  occurs  whenever  the  elements  Kmnm„  of  (10)  are 
always  near  unity,  as  they  are  for  practical  configurations 
producing  SBS.  To  see  this  we  consider  first  a class  of  pump 


waves  for  which  Eq.  ( 1 1 ) can  be  solved  exactly,  and  then  apply 
general  perturhation  theory  to  an  even  wider  class  of  pump 
waves. 

VI.  SET  OF  EXACT  SCATTERED  SOLUTIONS 
AND  THEIR  GAINS 

There  is  an  important  class  of  incident-mode  patterns  in 
a rectangular  waveguide  for  which  the  eigenvectors  of  both 
Hnm  and  are  tha  »ame  set,  provided  Kmnm„  — 1 las  for 
SBS)  which  we  assume  here.  This  is  th*  class  of  waves  for 
which  all  N excitable  modes  of  one  polarization  of  the  guide 
are  excited  with  “dual  energy  but  with  arhitrary  phases  <t>„, 

a„  -N- 1/2  *>-•>«.  (26) 

For  simplicity,  assume  the  first  M modes  in  both  x and  y 
directions  are  excited  (A/  * M-)  and  label  the  modes  by  x and 
y indices,  (n,,  ny ) *+  n such  that 

nx,  ny  * 0, 1,  2 M l (27) 

in  order  of  increasing  numbers  of  nodes.’  The  integrals  in  (15) 
and  (16)  are  then  easily  approximated  (assume  the  mode 
functions  are  sine  waves  vanishing  at  the  boundary)  to  ob- 
tain 

«../  - d„(  - [2  + Unx,lx  )|  1 2 + «(nv,/v)]/4S,  (28) 

where  the  4 function  is  1 if  its  arguments  are  equal,  and  zero 
otherwise.  It  is  easily  verified  by  substitution  in  Kq.  (11)  that 
the  N * M2  eigenvectors  have  components  (mode  ampli- 
tudes) 

Bn((*.  (>)  ■ M~x  exp  [i>„  + 2ri\nxlx  + nyly)/M\,  (29) 

where  pairs  of  the  integers 

(x.^-0,1 Ml  (30) 

are  convenient  labels  for  the  N vectors.  The  eigenvector  with 
l,  » ly  m 0 is  the  desired  phase-conjugate  state  (fi„  » a’„),  and 
it  has  the  largest  eigenvalue: 

Yoo  - a + 3.  (31) 

Then  there  are  2M-2  eigenvectors  having  either  l,  m 0 or  ly 
■ 0 (but  not  both)  which  all  have  the  eigenvalue 

Vo,  - a + (1  + 2Af)/4SM2.  (32) 

I 

Finally  the  (M-l)2  states  for  which  (r  ^ Oand  ly  9*  0 have  the 
eigenvalues 

y„„  - a + 1/4SM2  (33) 

For  the  sine-wave  modes, 

a ■ d * (1  + M~l  + M”2/4)/S.  (34) 

Therefore  the  “gap”  in  gain  between  Yoo  of  the  conjugate 
wave  and  Yo,  the  next  nearest  gain  is 

Yoo-  Yo.-(l-,&N-,/2)/S  (35) 

which,  for  a large  number  N of  modes,  approaches  the  gap 
between  Eqs,  (22)  and  (23)  for  Vm„  ■ 0.  Also,  for  large  i V,  we 
see  from  (21)  that  the  gain  of  the  conjugate  wave  is  the  same 
as  for  a plane-wave  pump  of  intensity  P/S  in  free  space. 

Studies  of  pump  beams  that  do  not  have  equal  power  in  all 
modes  have  led  us  to  believe  that  the  gap  between  the  highest 
and  next-highest  mode  gains  is  quits  generally  a less  a term 
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of  order  N~in  smaller.  This  gap  ia  usually  enough  in  practice 
to  make  a tingle  stimulated  solution  dominate  in  experiments 
where  the  backscattered  wave  grows  from  noise.  It  remains 
to  be  seen  in  Sec  VII  how  closely  this  tingle  solution  ap- 
proximates  a phase  conjugate  to  the  incident  wave. 


VII.  PERTURBATION  SOLUTIONS  FOR 
SCATTERING  FROM  ARBITRARY  PUMP  WAVES 

Having  seen  from  the  exact  solutions  for  B„  in  the  special 
case  of  the  previous  section  that  the  effect  of  V'mn  in  Eq.  (11 ) 
was  small  when  the  number  N of  pump-wave  modes  excited 
was  large,  we  are  encouraged  to  treat  Vmn  as  a small  pertur- 
bation on  Hmn  in  solving  the  eigenvalue  equation  (11)  (in 
which  we  continue  to  assume  Umn  - 0 for  reasons  argued 
previously).  Let  us  denote  the  exact  eigenvector  of  Hm„  + 
Vmn,  having  largest  eigenvalue  Vo,  by  Ron  Since  we  have 
already  shown  in  (2*2)  that  the  eigenvector  6J,01  of  Hmn  having 
highest  eigenvalue  C<01  - (a  -f  jj)  was  a*.  let  us  write  the  exact 
solution  as 

Bon-oi+c,.  (36) 

Then  standard  perturbation  theory  gives  for  the  correction 
c„  (repeated  mode  indices  are  to  be  summed  henceforth) 

Cn-  I 6Jr»&r  V,ma;/(G«»-G<»>)  (37) 

» *»  o 

to  lowest  order  in  the  perturbation  Vi„  Recall  that  the  bj,'1 
is  the  eigenvector  of  Hm„  with  eigenvalue  G*'1.  From  Eqs. 
(22)  and  (23)  we  recall  that  the  energy  denominator  G<01  - G'"1 
in  Eq.  (37)  is  a constant  0.  Since,  from  its  definition  in  Eqs. 
(13)-(17),  V„„  has  the  property 

amVmKa'n-0,  (38) 

the  sum  in  (37)  may  be  extended  over  all  » (with  constant 
denominator),  and  closure  invoked  to  obtain  the  expression 

c„“  (39) 

which  is  simpler  for  calculation.  The  most  important  quan- 
tity to  calculate  for  our  purposes  is  the  fraction  f of  power  in 
the  highest  gain  solution  Ron  that  is  not  in  the  phase-conjugate 
wave  a*.  Since  a\  a„  ■ 1 by  '!<), 

i - r/(l-r),  (40) 

where  r » cB  c„.  Fr  jm  Eqs.  (39)  and  (13)  we  have 

r-|o,|,«(Bl|am|*«mB|0n|Vd*  (41) 

where  ■ a,m  + 0tm  - o - 0.  This  is  the  most  important 
result  of  perturbation  theo.;*  for  assessing  the  fraction  f of 
stimulated  backward  scattering  .‘'at  is  not  phase  conjugate 
to  the  pump  wave.  If  r is  not  much  less  than  1,  the  process  is 
useless  and  it  is  pointless  to  calculate  further  corrections  to 
Eq.  (39).  Jf  r « l then  (39)  and  (40)  are  accurate  enough.  It 
is  of  minor  interest  to  note  that  the  gain  Y0  of  the  important 
wave  is  somewhat  larger  than  G,0>.  From  nondegenerate 
perturbation  theory  for  eigenvalues, 

Yo  - G«»  + rg  (42) 

to  second  order  in  Vmn.  (The  correction  linear, in  van- 
ishes because  of  Eq.  (38). ) That  is,  the  gain  is  increased  by 
the  fraction  r of  the  gain  “gap"  0. 


We  next  use  (41)  to  calculate  r for  a fairly  general  class  of 
putnp  beams. 

VIM.  MODEL  CALCULATION  OF 
NONCONJUGATED  POWER  FRACTION 

We  calculate  here  the  nonconjugated  power  ratio  r of  Eq. 
(41)  for  a realistic  set  of  complex  pump-beam  amplitudes  a*, 
each  member  having  arbitrary  phase  and  one  of  four  possible 
magnitudes,  in  the  ideal  rectangular  waveguide  considered 
in  Sec.  VI.  Expressing  the  mode  label  n in  terms  of  the  x and 
y mode-integers  l,  and  fv  of  Eq.  (30),  we  take  for  the  square 
magnitudes  required  in  Eq.  (41) 

kl*  -*,(/,)*,(/,).  (43) 

The  real  positive  functions  g,  and  gy  are  defined  to  have  the 
value  I*  for  s number  N+  of  the  Af  values  of  their  integer 
argument,  to  have  the  value  g.  (less  than  g+)  for  N-  of  these 
Af  values,  and  to  have  the  value  zero  for  the  remainder.  The 
particular  members  of  each  of  these  three  seta  need  not  be  the 
same  for  g,  and  gy.  This  results  in  four  possible  values  for 

kl*. 

g*S+.  g-g-,  t*g-  orO 

depending  orin  *•  The  total  number  N of  excitable 

modes  equals  Af1  (SN+N-). 

Again  we  assume  all  modes  to  have  the  same  linear  polar- 
ization („  parallel  to  t or  $ so  that  a„m  ■ 0„m  in  8nm  are  given 
by  Eq.  (28).  The  required  sums  in  (41)  are  then  easily  done 
to  yield 

r-4|(u*>*/(u*>«-ll.  (44) 

where  the  average  ( > is  to  be  performed  over  the  two-valued 
function 

+ + (45) 

with  the  two-valued,  normalized  probability 

Wi*g±NJ{N+g+  + N-g-).  (46) 

Clearly  r in  Eq.  (44)  is  a function  of  three  independent  pa- 
rameters, which  we  choose  as  N-,  and  g-/g+,  whose 
ranges  are 

0 S N*  S Af,  0 < N_  S Af,  OS  gjg+  S 1.  (47) 

Let  us  consider  r for  various  possible  parameter  seta  rep- 
resenting various  possible  classes  of  pump  beams. 

A.  N . (org~)  ■ 0;  1 < N+  S Af.  This  is  a pump  beam  with 
N * modes  excited  with  equal  amplitude  and  arbitrary  phase. 
In  this  case  Eq.  (44)  reduces  immediately  to  r ■ 0.  To  this 
order  in  perturbation  theory,  the  phase-conjugate  wave  with 
Bom  “ a*  is  exact.  We  saw  above  that  when  N+  “ Af,  this 
solution  is  exact  to  all  orders  in  V„m.  In  any  event  we  expect 
this  case  to  yield  maximum  gain  for  a wave  thet  is  indistin- 
guishable in  practice  from  the  phase  conjugate. 

B . g_*g+.  This  is  obviously  equivalent  to  case  A and  gives 

r ■ 0. 

C.  N-  ■ 0,  N+ ■ 1.  This  is  the  single-mode  pump,  which  we 
saw  above  had  the  conjugate  wave  as  an  exact  solution  to  Eq. 
(11);  and  Eq.  (47)  gives  r ■ 0 as  expected. 

D.  N«.*N-*1.  This  is  the  worst  case  in  the  model  (four 
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modes  with  three  amplitudes),  giving  the  largest  nonconjugate 
ratio  r for  given  g -/#♦  Simple  numerical  analysis  shows  that 
the  maximum  of  r is 

~ 0.0669  at  g . ~ 0.20g  ♦ . (48) 

That  is,  in  the  worst  case,  less  than  8%  of  the  hackscnttered 
energy  should  be  different  from  the  pure  phase-conjugate 

wave. 

E.  N+  " N-  » 1.  For  many  modes  excited  among  three 
amplitudes  this  way,  r N * a.  Plots  of  Eq.  (44)  show  that  the 
maximum  r is 

rm„~0.1N;2  (49) 

with  little  change  in  the"  range  0.1  < g-lg*  < 0.2  and  rapid 
drop  off  when  this  ratio  is  smaller  than  0.05  or  larger  than  0.5. 
A more  exact  treatment  is  not  useful  as  clearly  the  noncon- 
jugate fraction  f ~ r becomes  insignificant  when  the  number 
of  excited  modes  is  large. 

F.  N+  * N_,  ,V+  » 1.  Then  r <*  N71  as  in  Eq.  (49). 

In  conclusion,  we  see  that  for  any  pump  wave  that  is  dis- 
tributed among  guide  modes  with  arbitrary  phase  and  any 
distribution  among  four  amplitudes  obeying  (48),  the  re  is 
never  more  than  8%  nonconjugated  power  in  the  stiminrted 
wave  with  maximum  gain  and  generally  much  less.  V.  e ore 
led  to  conclude  (contrary  to  Ref.  4)  that  the  stimulated 
backscattering  wave  with  the  highest  gain  has  an  order-of- 
magnitude  more  power  in  the  phase  conjugate  than  in  the 
useless  background,  for  any  distribution  of  pump-mode  am- 
plitudes whatever,  provided  that,  as  in  SBS,  »>-u>  « w so  that 
Kmmrt  1-  We  consider  next  the  relaxation  of  this  condi- 
tion. 

IX.  PHASE-CONJUGATION  IN  STIMULATED 
RAMAN  SCATTERING  (SRS) 

When  the  scattered  frequency  w is  significantly  different 
from  »,  then  Kmnmn  in  Eq.  ( 16)  may  not  be  nearly  unity,  as  we 
have  assumed.  We  will  find  that  phase  conjugation  disap- 
pears as  |Km„mn  | declines,  but  this  need  not  happen,  even  at 
sizable  Stokes  shifts,  provided  that  the  guide  interaction 
length  L is  short  enough.  Before  deriving  the  conditions  w, 
»,  and  L must  satisfy,  we  consider  the  possibility  of  stimulated 
forward  scattering  (which  does  not  occur  in  SBS). 

Ifthefin  in  Eq.  (11 ) were  to  represent  amplitudes  of  forward 
traveling  waves  (i.e.,  -*  ~knj,  then  one  sees  immediately 

that  there  is  again  often  one  solution  having  about  twice  the 
gain  of  any  other  (when  Kmmn„  ~ 1),  but  that  solution  has 
~ o„.  That  is,  this  solution  is  a sort  of  replica  of  the  incident 
wave,  but  not  time  reversed  or  phase  conjugated.  This  pos- 
sibility may  be  of  interest  but  we  shall  not  consider  it  further 
here. 

To  consider  the  effect  of  reduced  Kmnmn  on  backward 
scattering,  we  consider  again  the  model  of  Sec.  VI  in  which  the 
first  y modes  along  both  the  x and  y axes  of  a rectangular 
guide  ore  excited  with  arbitrary  phases.  Inasmuch  as  this 
model  produced  results  close  to  that  of  the  variable-amplitude 
model  when  Kmnmn  ~ 1,  we  expect  it  to  do  so  here  also.  We 


specialize  to  a square  waveguide  for  which  the  transverse  ei- 
genvalues u„  of  Eq.  (7)  are 

u\  ■ *2(«;  + n;)/S  (50) 

with  the  mode  integers  as  in  Eq.  (27).  It  is  sufficiently  accu 
rate  to  expand  AA  of  Eq.  ( 10)  to  lowest  order  in  the  u • to  oh- 
tain  for  the  (m,n)  element 

AA  * q (mjf  -n;  4-  m*  - n;),  (51) 

where 

q - *AA/4 S (52) 

and  AX  is  the  difference  between  the  wavelengths  in  the  me- 
dium at  wand  *: 

AX  x (c/2x)  (1  /n„w  - l/n,i>).  (53) 

Here  and  n , are  the  refractive  indices  at  w anc  i\ 

We  use  these  specific  forms  now  to  calculate  the  noncon- 
jugated fraction  r,  using  the  perturbation  expression  (41),  and 
to  find  the  conditions  necessary  to  keep  this  he  low  some  de- 
sired maximum  value  tq.  With  Eqs.  (50)— (53)  in  (10)  and  ( 16), 
the  dmn  matrix  we  will  need  in  this  calculation  is  given  by 

4 Sdmn  >«  f dr[2e*«,,m$  “ "S'  + 4 (m,,  n,)| 

X |2e“»,,mJ  - "$»  + 6(mv,  ny ))//.,  (54) 

which  is  seen  to  reduce  to  the  previous  form  (28)  when  q — 
0.  The  <»ml  terms  in  (41)  for  r cancel  exactly  as  before,  leaving 
only  the  f)mn  terms  to  be  estimated.  This  we  do  by  calculating 
the  correction  A r that  is  lowest  order  in  qL,  by  expanding  (54) 
in  powers  of  q.  Higher-order  corrections  are  not  interesting 
as  they  only  tell  more  precisely  what  happens  after  the 
phase-conjugated  fraction  has  become  uninterestingly  small. 
We  will  also  approximate  the  sums  over  integers  in  the  re- 
sulting expressions  by  integrals,  since  the  accuracy  of  this  is 
quite  good,  especially  at  large  H . This  gives  directly  for  the 
increase  Ar  in  r from  its  value  when  K„nmn  ~ 1: 

Ar  » ^ (^f~L-)2(1  + 0 + 0 |M-'|).  (55) 

leading  to  the  condition  on  the  interaction  length  L 

L < 6 rln  S/N±\  (56) 

in  order  that  the  nonconjugated  fraction  be  kept  less  than  r0 
for  N equally  excited  modes  in  a rectangular  waveguide  of  area 
S.  We  feel  that  Eq.  (56)  can  be  applied  usefully  when  there 
are  N unequally  excited  modes  of  arbitrary  phase,  provided 
that  the  required  nonconjugated  fraction  r0  is  small. 

As  A r increases  by  (55),  so  does  the  gap  in  gain  between  the 
near'y-phase-conjugate  wave  and  its  nearest  competitor  de- 
crease. An  estimate  of  this  can  be  made  from  Eqs.  (22)  and 
(23).  However,  it  is  not  very  interesting  to  calculate  this  gap 
when  the  waves  are  already  degraded. 

In  the  SRS  phase-conjugation  experiment  of  Zel'dovich  et 
oi., s the  incident  beam  was  focused  into  an  "infinite"  medium 
of  CS^.  This  was  roughly  equivalent  to  employing  a wave- 
guide whose  length  was  a diffraction  length  -'SX_l  N~w2  a 
We  may  use  Eq.  (55)  then  to  estimate  the  nonconjugate  frac- 
tion A r arising  from  this  finite  length,  obtaining  Ar  ~ 0.7"i>, 
in  agreement  with  the  qualitative  report  of  excellent  conju- 
gation. (See  parameter  summary  above.) 
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In  conclusion,  stimulated  backward  scattering  at  u in  a 
waveguide  can  be  an  efficient  generator  of  the  phase  conjugate 
of  an  incident  wave  at  r,  irregardless  of  the  distribution  of  the 
N mode  amplitudes  or  the  incident  wave,  provided  that  the 
frequency  u>  it  not  so  differ  ::n  from  r for  given  guide  length 
as  to  violate  Eq.  (56).  For  some  incident  beams  the  noncon- 
jugate fraction  contained  in  the  backward  wave  with  highest 
stimulated  gain  may  be  of  order  10-1  but  is  generally  much 
smaller,  decreasing  as  N'1  as  the  number  of  modes  N in- 
creases. Scattering  without  a waveguide  may  be  approxi- 
mated by  considering  a waveguide  one  diffraction  length  long. 
However,  when  a waveguide  will  satisfy  the  above  criteria  for 
phase  conjugation,  it  has  the  advantage  of  greatly  reducing 
power  requirements  and  the  competition  from  other  nonlinear 
effects,  such  as  self- focusing. 
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7The  degeneracies  referred  to  here  are  those  arising  from  symmetry, 
such  as  occur  for  similar  right-  and  left-circularly  polarized  modes 
in  a cylindrical  waveguide.  That  AA  cannot  approach  zero  (to 
within  less  than  L~l > for  any  other  cases  is  because  otherwise  u, 
4 u}  cannot  come  within  S' 1 of  uj|  + uf  without,  at  the  same  time, 
to, 're  being  a drastic  reduction  in  the  z-y  integral  in  (9).  This  may 
be  appreciated  by  studying  the  example  of  Eq.  (50)  for  uj.  Here, 
without  the  aforementioned  degeneracies  occuring,  one  sees  that 
AAL  i L/Sk,  and  L/Sk  is  never  much  less  than  unity  when  guiding 
occurs. 

*A  waveguide  which  imitates  well  the  interaction  of  focused  unguided 
waves  may  be  imagined  as  follows.  Construct  that  complete  or- 
thonormal set  of  free-space  Guassian-beam  modes  whose  param- 
eters are  such  that  the  smallert  number  N of  modes  need  be  su- 
perposed to  give  a good  representation  of  the  (focused)  incident 
bea'i.  Then  let  the  waveguide  axis  coincide  with  the  z axis  of  these 
modes  and  let  it  barely  encompass  the  beam  waist  uver  the  length 
where  the  waist  size  does  not  change  appreciably.  This  length  it 
generally  of  order  SA  “ 1 N * ,/J.  Calculate  the  backscattered  wave 
using  the  mode  decompositions  in  Eqs.  (5)  and  (8)  is  we  have  pre- 
scribed. The  field  patterns  calculated  at  the  entrance  to  the 
waveguide  should  approximate  those  in  the  same  plane  in  free 
space,  since  guiding  is  minimal. 
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ABSTRACT 


In  current  experiments  studying  CW  optical  wavefront  conjugation  by  degenerate 
four  wave  mixing,  the  of  foots  of  atomic  motion  are  not  negligible.  1 summarize 
a calculation,  that  includes  such  effects  from  the  beginning,  of  the  small 
signal  phase  conjugate  reflection  coefficient  for  both  Doppler  broadened  and 
homogeneously  broadened  resonant  transitions. 
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Wave  Propagation  Lab 
Boulder,  Colorado  80505 
(505)  499-1000 


111 


- 

PRECEDING  PaQ£  &ULMK 


In  wavefront  conjugation  experiments  employ  inn  degenerate  four  wave  mixing 
in  an  atomic  vapor,  atomic  motion  can  have  a non-negl igiblc  effect  even  in  a 
homogeneously  broadened  medium  if  the  longitudinal  (population  difference) 
relaxation  time,  T.  , is  longer  than  the  characteristic  Doppler  time  (optical 
wavelength/thermal  speed).  The  form  of  the  solution  given  in  the  analyses  of 
previous  authors*’2  depends  only  on  the  homogeneity  and  isotropy  of  the  medium  and 
the  validity  of  the  KKB  approximation  (slowly  varying  amplitudes);  therefore  the 
effects  of  the  atomic  motion  arc  totally  contained  in  the  coupling  , ammeters 
ol  and  ts , that  respectively  express  t lie  polarizability  and  the  coupling  of  the 
forward  and  backward  waves,  in  the  presence  of  counterpropagat ing  pump  beams. 
Because  the  phase  conjugate  reflection  coefficient  R is  proportional  to  the 
fourth  power  of  the  pump  field2,  a lowest  (third)  order  perturbation  theory 
calculation  need  only  consider  the  coupling  constant  k.  In  this  order  we  thus 
have 

R = — [1  -exp(-2a  L)]fc  , U) 

4a2 

r 

where  a is  the  linear  absorption  coefficient  and  L is  the  length  of  t he  inter- 
r 

action  region  along  the  direction  of  the  signal  beam. 

In  order  to  calculate  k,  including  the  effects  of  atomic  motion,  go  back 
to  the  Bloch  equations3  ’ 4 for  the  two -by -two  dd.  .'ty  matrix,  written  "symbolically" 
as 

Iy  P(t)  = - £ [HCt) .Pit)]  (2) 

+ relaxation  (Ti.Tz)  terms 


where 


hdln 


II  = 


2 ^0 


■ d,  , ,0  I-. 

) - pb (t.)  (1  w) 


(3) 


112 


(Wo  is  the  resonant  frequency  of  the  transition,  |i  the  dipole  matrix  element 
and  li  the  electric  field.)  The  position  and  \clocity  dependence  of  the  density 
matrix  can  then  be  taken  into  account  by  making  the  replacements 


P(t)  + p(x,v,t) 
E(t)  + E(x,t) 


9_ 

9t 


d_ 

dt 


9 

9t 


+ 


. 


I have  found  it  convenient  to  employ  the  vector  representation  of  Feynman,  Vernon 
and  Hcllwarth,5  in  whic!  the  resultant  equations  are 

( ft  + ^ + { P 3 (x , v , t ) - p3(o)(v)}  = ^-E(x,t)p2  (4) 

( ft  + ^ + f7  ^ P»  Cx,  v, t)  - O)0  P2  (5) 

( It  + ^ + T7  ' P?i>  'V,t^  = 100  Pl  + 17"  E("^» fc)P3  (6) 


I take 


Po 


(71  U2)3/2 


exp( ) 

u2 


(7) 


where  u is  the  thermal  speed. 

These  eq"- tions  can  he  derived  somewhat  more  rigorously  by  considering  the 

"quasi-classi  il"  (•&>())  limit  of  the  equation  obeyed  by  the  Wigner  function.  An 

additional  benefit  of  the  latter  approach  is  the  presence  of  terms  which  describe 

velocity  changing  (recoil)  effects,  corresponding  to  the  F • — terms  in  the 

9v 

classical  Boltzmann  equation.  The  inclusion  of  these  terms  will  necessitate  the 
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introduction  of  yet  another  phenomenological  relaxation  time  for  the  trace 
of  the  density  matrix  (number  density).6  The  full  equations  could  be  used 
to  quantitatively  explore  the  recently  proposed7  idea  to  use  density  changing 
effects  to  enhance  nonlinear  couplings  in  at'  ••  vapors. 

The  polarization  is  given  by 

P(x,t)  = p / d3  v pi(x,v,t)  . 


The  third  order  perturbation  theory  result  for  pi  is  most  easily  obtained  by 
first  going  to  an  interaction  representation  (in  which  vanishes  in  the 

absence  of  an  electric  field),  Fourier  transforming  in  both  space  and  time 
(giving  integral  equations)  and  finally  dropping  non-resonant  terms  (the 
"rotating  wave"  approximation).  In  order  to  perform  the  velocity  integration 
in  the  above  expression  for  p(x,t),  Feynman's  trick8  of  combining  the 
perturbation  theory  denominators  is  invaluable.  The  phase  conjugate  coupling 
amplitude  k is  then  identified  by  picking  out  the  part  of  P (q)  proportional 
to  l:*(q),  where  q is  the  wave  vector  of  the  signal  wave.  The  result  is 


Y1Y2  l'  1 

K*  = do J dx  J dy 


n0 

I-2 

s 

1-x 


f 3 (x,y,0) 


0 0 

((1-x)  (i  + 6)Y2  + ixYi 
f(x,y,9) 


ri-1--  1 j...  7i  / 

((1-x) (i*6)Y2+ixYi  ^ 

Li  6>'2  d2  (x  ,0)  ' 

k d (x,0)  • / 

| ♦ {Q  - 0}] 


(9) 


where 


Yi  ,2 


J 1 
T 1 ,'2  ku 


(10) 
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d(x,0)  = / 1-2XCOS0+X7 


(ID 


t'U.y.O)  - (l-2y)zcosJ~  + (l+x)?sin?  | > (12) 

0 is  thy  angle  between  the  signal  and  the  pump  beams,  6 is  the  normalized 
(to  1/T2)  detuning  from  icsonancc,  and  Z'  and  Z"  are  the  first  and  second 
derivatives  of  the  plasma  dispersion  I unction. 9 In  order  to  proceed 
further,  approximations  to  Lite  remaining  (leyniimu  parameter)  integrations  are 
used  appropriate  to  either  the  homogeneously  broadened  regime  ( y ? > > I ) or  the 
Doppler  broadened  regime  (y2<<l). 

In  the  homogeneously -broadened  regime,  the  important  region  ot  integration 
is  x~l.  Approximating 

d ( x , 0 ) :i  d(l  ,0)  (13) 

f(x.y.G)  = f ( 1 , y ,0)  (14) 

and  integrating  by  parts,  one  obtains 

2 n0  I i i 1 i 2 

k = — — mfn.O)  , (IS) 

I 1 1 m " ) ( i i A ) I.'’ 


where 


m(Yl,0)  - ? (s  + cxp(s+)erfc(sj)  + s_exp(s2_)  crfc(s_)J 


(1(») 
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(17) 


Y 

i 

2 sin  0/2 
Y) 

2 tos  0/2  (18) 


Ihus  the  k of  reference  1 is  multiplied  by  the  attenuation  factor  in,  and  K hy 
its  square.  Iliis  lactor  is  plotted  in  figure  1 as  a function  of  angle  for  constant 
Y i The  asymtotic  value  m(0=0)  -►  'j  for  infinite  thermal  speed  (Yi  0)  is  a 
reflection  of  the  fact  that  only  one  of  the  two  possible  holographic  gratings 
formed  by  the  signal  and  one  pump  beam  (the  infinite  wavelength  grating)  survives 
the  "washing  out"  by  the  atomic  motion. 

In  the  Doppler  broadened  regime  (Y2«l),  the  derivatives  of  the  Z function 

can  be  replaced  by  their  small  argument  limits,  as  long  as  0 does  not  become  so 

small  as  to  makdc  Yi,2/f  or  Yi.2/d  large.  The  result  under  these  circumstances 

is  2 

liji-2  V2 

k s 2i  /n  a0  — ; 0 / 0 

hji  (ID) 


(Ibis  is  only  valid  if  the  detuning  uoes  not  exceed  t he  Doppler  line  width, 
6y2«1.)  In  order  to  estimate  the  behavior  for  small  angles,  < can  be  cal- 
culated for  O = 0: 


K 


% jJ'j>  Y2 

2 bs2  i+6 


0 = 0 


(20) 


comparing  these  two,  we  see  that  the  small  angle  behavior  sets  !■  it  an  angle 
0q , given  by 
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sin1*  Oo 


2 > 

10  Y , Y ■/  (1*.S2) 


(21) 


For  fixed  G,  the  reflection  coefficient  will  exhibit  a linewidth 


sin"  n 

(A'o)2  = min(rd2,  r Pd2) 


(22) 


where  F , ku.  A simple  interpolating  fornnila  for  the  Doppler-broadened 
reflection  coefficient  that  agrees  w.f  n the  above  small  and  larg^  angle 
results  is 


2 2 

Y Y 

1 2 


II  2 


sin  0 + lo 


2 2 

Y Y 
1 2 


(1  + 62) V H 2 / 


1 1 -exp( -2«rL) J ; 
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This  is  plotted  (or  fixed  detuning  parameter  6 in  figure  2 and  for  fixed  angle 
in  Figure  3.  The  angular  dependence  is  qualitatively  similar  to  the  homogeneously 
broadened  case,  the  main  new  feature  arising  being  the  quasi  l)oppler-free  nature 
of  the  frequency  dependence  at  small  angles,  consistent  with  the  observations 
of  Liao,  Bloom  and  liconomou . 1 0 
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Figure  Captions 


Figure  1 


Figure  2 


Figure  3. 


Amplitude  reduction  factor  m(Yi,0)  versus  angle 
broadened  system) . 


(for  homogeneously 


Interpolated  reflectivity  versus  angle  (normalized  to  0*6=0)  for 
Doppler  broadened  system.  Yi  y2  = jo'* 

Reflectivity  lineshnpes  (Doppler  broadened  system)  for  various  angles 
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Figure  3.  Reflectivity  lineshapes  (doppler  broadened 
system)  for  various  angles. 
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APPENDIX  D 


THE  IMPACT  OF  IS0PLANAT1C  EFFECTS  ON  TARGET  REFERENC  INC 
SYSTEMS  AT  VISIBLE  WAVELENGTH* 

In  this  class  of  problem  the  transmitted  laser  radiation  Is  to  be 
focused  to  a spot  on  the  target.  The  reference  information  is  from 
laser  reflection  off  of  a glint  or  highlight  located  at  or  near  this 
spot.  Specifically  with  orbiting  targets,  as  illustrated  in  Figure 
D-l,  the  target  emits  or  reflects  a reference  wave  at  position  A on  its 
orbital  path  and  then  moves  to  position  B during  the  round  trip  propa- 
gation time.  That  is,  it  is  required  to  point  the  transmitted  laser 
beam  along  a path  which  intercepts  the  target  at  a point  slightly  ahead 
of  its  position  at  the  time  the  laser  energy  leaves  the  transmitter 
aperture.  Thus  the  laser  must  be  fired  ahead  of  the  target  in  order 
to  intercept  it.  For  this  application,  the  reference  beam  and  the 
laser  beam  are  essentially  similar  cones  with  axes  inclined  at  the 
point-ahead  angle  defined  by  the  angular  distance  between  the  positions 
A and  B.  For  low  orbits  the  xequired  point-ahead  angle  is 

2v0 

0 ■ cos  , (D.  1) 

pc 

where  v is  the  orbital  velocity  (typically  of  the  order  7 km/sec  for 
low  altitude  satellites),  c is  the  velocity  of  light  and  $ is  the 
azimuthal  angle.  At  zenith,  a typical  point-ahead  angle  is 


0 % % 50  prad. 

P 3x10® 


CD.  2 > 


The  work  described  in  this  appendix  was  performed  under  IR&D  and  not  in 
the  course  of  the  contract.  It  is  included  because  of  its  relevance 
to  atmospheric  phase  compensation  with  nonlinear  phase  conjugation. 
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In  spite  of  the  fact  that  a synchronous  satellite  remains  essentially 
fixed  at  a given  point  above  the  surface  of  the  earth,  there  is  still 
a point-ahead  requirement  in  this  case  also.  The  earth  and  the  satel- 
lite rotate  in  unison  at  the  rate  of  2tt  radians  per  24  hours.  Thus, 
the  reference  (transmitter)  aperture  must  be  pointed  at  the  position 
where  the  transmitter  (target)  aperture  will  have  moved  to  during  the 
time  it  takes  the  light  to  propagate  through  the  intervening  space. 

For  synchronous  orbits,  the  relation  in  Eq.  D.l  remains  valid  if  we 
replace  the  cosine  factor  by  unity  and  use  the  value  SIR  for  the  orbital 
velicity  where  Q is  the  rotation  of  the  earth  and  R is  the  altitude  of 
the  orbit.  This  yields 


» 2X(2tt/3600)  X 4 X 107 
3 X 108 

* 20  yrad.  (D.3) 

The  conventional  isoplanatic  problem  in  which  the  reference  and 
transmitter  beams  have  similar  geometries  has  been  widely  studied  with 
two  different  approaches.  The  problem  was  first  studied  by  Fried 
using  a Huygens-Fresnel  type  of  calculation.  Fried's  theory  requires 
the  evaluation  of  some  fairly  complicated  integrals  which  have  been 
coded  for  numerical  evaluation  by  personnel  at  the  Air  Force  Weapons 
laboratory  (AFWL).  Using  the  AFWL  code,  we  have  recently  obtained 
results  from  the  Fried  theory  for  isoplanatic  effects  at  3.8  microns 
as  a part  of  a ground-based  beam  control  study  presently  being  perform- 
ed for  the  Air  Force.  In  addition,  we  have  modified  our  propagation 
and  adaptive  optics  simulation  software  to  allow  us  to  model  the 
Isoplanatic  problem  without  making  the  simplifying  assumptions  inherent 
in  the  Fried  theory.  In  each  of  these  approaches  we  have  assumed  an 
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index  structure  constant  variation  with  altitude  with  a similar  to  the 
model  proposed  by  Hufnagel7  in  1974.  It  consists  of  the  1 974  liufnagel 
model  plus  an  additional  ground  layer  term  which  is  used  to  adjust  the 
strength  of  this  layer.  Hufnagel's  model  has  the  form 


CN2  - 2.  7[2. 2(  102*3h) 10  expf-hX^yr 


10'16  exp(-h/ 1.5)),  m_2/3 


( D . A ) 


where  h is  the  altitude  in  km  and  v is  the  rms  wind  speed  between  5 

w * 

and  20  km.  In  our  studies  we  have  used  three  different  values  of  v 

w 

(15,27  and  40  m/sec). 

Figure  3-2  shows  a comparison  of  results  obtained  from  the  Hughes 

propagation  code  and  from  the  Fried  theory.  These  results  apply  to 

propagation  at  3.8  microns  and  assume  perfect  adaptive  optics;  i.e., 

infinite  temporal  and  spatial  bandwidth.  In  these  calculations  we 

used  the  Hufnage 1 model  given  in  Eq.  D.4  (v  - 40  m/sec)  augmented  by 

w 

a ground  layer  of  the  form 

<Cn3)  8™md  ' ‘O'14  «*P<-10h),  i"'2'3  ■ <».5> 

The  solid  curve  in  ^igure  D-2  is  the  Fried  theory  result.  The  circles 
indicate  the  average  result  obtained  from  the  propagation  code  and  the 
bars  indicate  the  variation  observed  about  this  average  (recall,  that 
the  propagation  code  generates  results  for  typical  sample  media  from 
the  ensemble  of  possible  media  and  thus  yields  results  which  inherently 
fluctuate).  The  agreement  between  the  results  in  this  case  is  strik- 
ingly good,  which  implies  that  the  assumptions  inherent  in  the  Field 
theory  are  valid  under  these  conditions.  Moreover,  these  results 
indicate  that  the  isoplanatic  effects  at  3.8  microns,  although  bother- 
some, can  probably  be  tolerated  for  the  typical  point-ahead  angles 
cited  in  Eqs.  D.2  and  D.3. 
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In  contrast  to  the  relatively  benign  nature  of  the  isoplanatic 

problem  at  3.8  microns,  the  decrease  in  Strehl  ratio  caused  by  the 

point-ahead  angles  given  in  Eqs.  D.2  and  D.3  is  quite  significant  at 

visible  wavelengths.  An  indication  of  the  severity  of  the  isoplanatic 

problem  at  .5  microns  is  given  in  Figure  D-3.  These  results  were 

obtained  from  results  given  by  Fried  in  his  1977  report.  They  are 

replotted  from  the  data  given  in  Figures  3a  and  3b  of  that  report  and 

apply  to  two  different  structure  constant  profiles.  One  of  these 

profiles  is  based  on  the  data  of  Miller  and  Zieske1  and  the  other  on 

2 

the  data  of  Barletti.  As  shown  in  Figure  D-4 , these  two  models  are 
similar  below  an  altitude  of  5 km  but  differ  significantly  above  that 
altitude.  The  Barletti  model  fits  the  Hufnagel  model  used  in  our 
3.8  micron  work  reasonably  well  for  a rms  wind  velocity  of  27  m/sec. 

We  believe  that  the  Barletti  data  are  probably  more  representative 
of  the  type  of  turbulence  profile  that  will  be  encountered  but  this 
question  must  be  explored  experimentally  to  verify  this  conslusion. 

For  either  profile,  however,  the  key  conclusion  is  thr.t  the  isoplanatic 
effects  range  from  severe  to  overwhelming.  In  order  to  obtain  a 
compensated  Strehl  ratio  greater  than  .8,  it  will  be  necessary  to  limit 
the  angle  between  the  reference  and  transmitter  beams  to  less  than  2 
microradians,  if  the  3arletti  data  apply,  and  to  less  than  6 microradians 
if  the  more  benign  data  of  Miller  and  Zieske  apply.  To  our  knowledge, 
the  only  way  to  achieve  this  will  be  to  use  a reference  which  is  locat- 
ed on  a "stick"  or  on  another  satellite,  which  leads  the  relay  satellite 
in  orbit. 

We  would  like  to  emphasize  that  the  above  conclusions  relative  to 
the  magnitude  of  the  isoplanatic  effects  at  visible  wavelengths  are 
based  entirely  on  results  obtained  from  Fried's  theory.  We  have  not 
used  the  Hughes  adaptive  optics  simulation  to  check  these  conclusions. 
Unfortunately,  this  software  must  be  modified  before  we  can  obtain 
results  pertinent  to  the  extremely  large  transmitter  apertures  contem- 
plated in  the  systems  of  interest. 
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Figure  D-l.  The  basic  isoplanatic  problem.  The  target  emits 
or  reflects  a reference  wave  (dashed)  at  position 
A and  moves  ahead  to  position  B during  a round- 
trip  propagation  time. 
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Figure  D-2.  The  conventional  isoplanatic  degradation 
experienced  at  3.8  pm  with  a 2.5  meter 
aperture.  Comparison  results  obtained 
using  two  basically  different  classes  of 
calculation  are  used  to  compute  the  strehl 
ratio  associated  with  a given  offset  angle. 
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Figure  D-3.  The  conventional  isoplanatic  degradation  for 

^ " . 5 ym  for  a 4.0  meter  aperture  after  Fried 
with  ideal  adaptive  optics. 
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ALTITUDE,  KM 

Figure  D,-4.  Measurement  of  the  vertical  distribution  of  the 
optical  strength  of  turbulence.  The  Barletti 
model  was  obtained  by  thermosonde  baloon.  The 
Miller  & Zieske  data  only  apply  above  4 km. 


130 


appendix  e 


DOPPLER  AND  POINT-AHEAD  OVERRIDE  TECHNIQUES  IN 
FOUR-WAVE-MIXER  PHASE  CONJUGATION 

For  an  orbiting  retroreflector  with  orbital  velocity  v and 
line-of-sight  angle*  0,  as  illustrated  in  Figure  E-l (a) , the  required 
point-ahead  angle,  as  seen  at  the  conjugation  subsystem,  is+ 


a 

pa 


cos 


e 


(E-l) 


whre  M is  the  optical 
typically  fall  in  the 
is 


magnification**  of  the  beam  director  (M  will 
range  40  < M < 400).  The  fractional  Doppler 


shif  t 


5 


I) 


2v 

— sin  0 

c 


(E-2) 


uhero  v Is  the  orbital  tangential  velocity,  and  C is  tile  velocity  of 
ligtt.  Clearly  „pa  and  differ  only  in  scale  (the  M multiplier)  and 
in  their  8 dependence.  The  rav  (M  - 1 ) point-ahead  angle  typically  is 
in  the  renge  of  40  to  60  grad. 

The  normal  k-vector  phase-matching  relations  are  given  in  Figure 
( ),  while  the  k vector  diagram  for  an  equal  frequency  (degenerate) 
four-wave  phase  conjugator  is  illustrated  in  Figure  E-l(c).  Since  the 
system  is  degenerate,  the  angle  0 between  pumps  and  signal  is  arbitrary. 


+F°r  low  orbits,  0 is  approximately  the  zenith  angle. 

The  sign  ambiguity  reflects  the  possibility  that  there  may  be  one  or 

r ?:„‘r'rs,ioi; in  the  °pticai  * ■»„  mdSe . 

the  inpit  beam  P°Sitl°"  ot  the  “"I-*1'1-"  relative  to’ 


** 


For  the  "carrot-on-the-stick" 
compensation  is  not  employed, 
compensation  will  be  absorbed 


approach,  where  full  point-ahead 
the  residual  point-ahead/-behind 
within  the  M factor. 
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<0  NORMAL  (DEGENERATE)  PHASE  CONJUGATION 

* VECTOR  DIAGRAM.  THE  ANGLE  p IS  ARBITRARV  AND 


However,  if  the  output  must  be  shifted  (slightly)  in  frequency  to 
compensate  ^ above,  the  degeneracy  is  formally  broken,  and  the  angle  U 

is  not  arbitrary. 

Figure  E-l(d)  illustrates  a k-matched  system  in  which  the  frequency 
of  pump  A?  is  shifted  upward  by  a small  fraction  fip2.  causing  a 
corresponding  fractional  shift  in  the  conjugated  output  A3  (the  original 
frequencies  are  normalized  to  unity  and  the  sizes  of  6^  and  ^ arc 
greatly  exaggerated  in  the  figure).  For  small  6 and  *.  the  geometry 

gives 

h - 1(1  A tan -§-  . (E-d) 

To  compensate  for  Doppler  and  point-ahead  simultaneously,  we  can  let 


Dividing  Eq. 
E-l  gives 


ijj  ■ lot 


pa 


6p2  " 6D 


(E-4a) 
(E-4b ) 


E-4a  by  Eq 


E-4b  and  substituting  from  Eqs.  E-3,  E-2,  and 


tan  j - ±M  cot  0 


(E-5n) 


or 


6 


it  + 2 tan 


(E- 5b) 


Typically,  operating 
for  typical  M values 
by 


ranges  of  0 will  fall  in  the  range  -n/4  < 0 

(40  < M <400),  Eq.  E-5b  can  be  well  approximated 


B 


TT 


¥ 


20 

M 


(E-5c) 
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With  M ■ 1,  even  a modest  range  of  (-45  < ft  < 45°)  requires  a wide 
swing  in  pump  orientation,  and  the  pumps  become  nearly  coi linear  with 
the  input-output  waves  at  0 - 0 causing  a separability  problem.  However, 
with  modest  magnifications  (M  ■■  40),  the  required  swing  in  pump  orienta- 
tion is  greatly  reduced,  although  the  separability  problem  remains. 

At  the  overhead  condition  (0  - 90°),  point  ahead  is  required  but 
not  Doppler  compensation.  However,  a second-order  fractional  frequency 
shift  6^2  i*'  PumP  is  required  to  produce  the  desired  point-ahead 

angle  <l> . This  shift  is  given  to  an  excellent  approximation  (with  an 
exact  k vector  match)  as 


3 

For  v - 7 x 10  km/sec,  M ■■  40,  and  0-0: 


p2 


. 4> 

\ 3 x 108  / 


,4  x 10 


(E-7) 


At  a wavelength  of  3 pm,  for  example,  this  corresponds  to  only  a modest 
frequency  shift  in  pump  2 of  about  44  MHz.  Another  possible  approach 
for  this  overhead  regime  (which  avoids  the  frequency  shift  requirement* 
entirely)  is  to  let 


k2  - k4  (E-8a) 

kj  - k'3  , (E-8b) 


as  illustrated  in  Figure  E-2. 


Actually,  this  is  only  a formal  requirement.  Practically  speaking, 
the  k vector  mismatch  resulting  from  no  frequency  shift  in  tolerable. 
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Up  to  this  point,  we  have  assumed  that  the  probe  wavefront  is 
a plane  wave.  Of  course,  for  applications  of  real  interest,  it  is  not. 
How  does  the  tilt/frequency  offset  affect  the  phase  conjugation  process? 
As  shown  below,  we  effectively  repoint  the  conjugated  wavefront  much  as 
a tilt  mirror  would  do. 

For  purposes  of  analysis,  we  have  defined  the  somewhat  simplified 
system  illustrated  in  Figure  E-3  in  which  the  return  wave  front  is 
essentially  propagating  in  the  negative  z direction  and  has  picked  up  a 
phase  distortion  4> (x ,y ) from  a phase  screen  reasonably  near  the  primary 
optics.  Thus,  we  assume  the  incident  field  through  the  collecting 
aperture  to  be* 


i<|>(x,y)  i(-k  z-w,t) 

* A^F(x.y)  e e , (E-9) 


where  F(x,y)  is  the  collecting  aperture.  After  demagnification,  the 
incident  field  on  the  conjugator  is  approximately 


i$(Mx,My)  i(-k  z-w  t) 

U7  1 ■ A,  1 e e » (E-lOa) 


where 


A.'  ■ MF(Mx,My)A.  (E-lOb) 

and  where  we  have  assumed  that  the  propagation-diffraction  effects 
through  the  telescope  are  minor.  The  pump  waves  are  assumed  to  be 
plane  waves: 


*It  will  be  understood  that  the  real  fields  are  the  real  parts  of 
the  given  expressions. 
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Figure  E-2.  Point-ahead  and  k vector  matching  without  a requirement 
for  pump  frequency  perturbations. 
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Figure  E-3.  Geometry  and  coordinate  system  for  analysis 
of  the  perturbed  system. 


(6-1 1 ) 


i(k.  -r-ui.t) 

Ul-Ale 

i(k  *r-w  t) 

U2  - A2  e J (E— 12) 

The  nonlinear  outputs  of  present  interest  result  from  the  third  order 
expansion  in  the  total  field: 

U - C (U.  + U0  + U.  + c.c.)3  , (E— 13) 

O 1 L 4 

where  c.c.  represents  the  complex  conjugate,  and  C is  a conversion 
efficiency  parameter.  In  particular,  the  conjugated  output  of  inter- 
est will  be  a component  called  U^': 

u3’  - c[u4  * u2u1  + c.c.] 

= CAjA^A^'  exp  (-!4>(Mx,My)  ] exp  i ( • r - wet)  , (E-Ha) 


where 


k3  * kj  + k2  - k4 


(E-Hb) 


and 


k,  * k - a (E-Hb) 

4 z z 

;o„  - w.  + to.  - w,  (E-Hc) 

3 2 14 
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and  where  »z  la  a unlt^  vector  ln  the  z direction.  In  the  offset  desl8n. 
w«  have  arranged  that  be  offset  from  by  , - * auch  that 


<3  • r - ZKz  (1  - 6)  cos  + yK2  (1  - 6)  sin  ^ . (E_15 

Further,  It  has  been  proposed  that  the  frequency  of  pump  one  equal  that 

of  the  return  wavefront  and  that  the  frequency  of  pump  two  be  perturbed 
from  the  return  wavefront  such  that 


“l  " “4  (E-I5b) 

and 


“2  " “4  O “ 6) 


(E-15c) 


giving 


“3  " (1  “ 6) 


(E— 1 5d ) 


Under  these  circumstances,  it  follows  that 


U3-  CA1A2A4  exP  [^(Mx.My)]  exp  i^cos  *+  y^.in  * - V ] (1-6) . (D-16) 

This  Is,  to  a first  approximation,  a pseudo  plane  wavefront  of  frequency 

propa8aiin8  at  a 8ma11  offset  angle  <p  from  the  z axis  and  having 
superimposed  on  it  the  demagnified  wavefront  error.  We  assume  that  the 
reverse  propagation  distances  to  the  primary  are  relatively  short  such 
that  his  field  propagates  like  a plane  wave  to  the  main  aperture,  retain 
ing  its  phase  distortions.  However,  the  demagnifications  are  stripped 
off,  giving  an  aperture  field  at  : • z which  is 

d 


Uj-CjA.A^Ftx.y)  e*p  C- 1* (x .y > ] exp  amp.  v,(1_s)  , (E„17a) 
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where  a,  a constant  of  no  importance  here,  is 

*a  " Dkz  cos  ♦ • (E-l/b) 

The  field  is  effectively  an  aperture  excitation  field  which 
reradiates  at  the  off-axis  angle 


a ■ ct 


pa 


sin  ijj  . 

M M ‘ 


>' E -18) 


However,  this  tilted  wavefront  reradiates  with  the  wavefront  error 

distribution  ~4> (x ,y ) superimposed  on  it.  Functionally,  the  resulting 

radiating  source  could  be  replaced  with  a thin  phase  screen  [ -<p (x , y ) ] 

in  front  of  the  aperture.  Small  changes  in  the  exit  angle  by  which  the 

exit  wave  passes  through  this  screen  do  not  change  its  effect  on  the 

radiating  field.  Thus,  for  small  angles  < p,  the  radiated  fields  are 

equivalent  to  what  the  reradiated  (conjugated)  wavefront  would  have 

been  without  the  pump  perturbations,  except  that  the  entire  radiation 

pattern  is  tiled  by  an  angle  a , and  the  reradiated  frequency  ui  is 

pa  3 4 

changed  with  a fractional  decrease  6. 

Of  course,  if  the  offset  angle  is  appreciable,  an  isoplanatic 
problem  exists,  as  discussed  in  Appendix  D.  If  the  frequency  offset  is 
too  great,  large  optical  path  differences  may  not  be  well  compensated. 
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